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PREFACE 

This report is one of seven separate reports prepared 
by six discipline-oriented analysis teams of the Earth 
Observations Division at the NASA Lyndon B. Johnson Space 
Center, Houston, Texas. 

The seven reports were prepared originally for Goddard 
Space Flight Center in con^liance with requironents for the 
Earth Resources Technology Satellite (ERTS-1) Investigation 
(ER600) . The project was approved and funded by NASA 
Headquarters in July 1972. 

This report (Volume IV) was con^iled by the Range 

Analysis Team. The following is a list of the teeua meiabers. 

E. P. Kan, Lockheed Electronics Company, Inc. 

J. M. Disler, Lockheed Electronics Company, Inc. 

H. D. Parker, Jr., Lockheed Electronics Company, Inc. 

A. Simon, Lockheed Electronics Company, Inc. 

D. Pendleton, Soil Conservation Service, USDA. 

The total investigation is documented in the following 
reports. 
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THE ERTS-1 INVESTIGATION (ER-600) 

VOLUME IV — ERTS-1 RANGE ANALYSIS 
(REPORT FOR PERIOD JULY 1972 - JUNE 1973) 

By R. Bryan Erb 
Lyndon B. Johnson Space Center 

1.0 SC^ MARY 

The Earth Resources Technology Satellite (ERTS-1 ) Range 
Analysis Teeun was formed in June 1972 to conduct a 1-year 
investigation of the utility of ERTS-1 remote sensor data in 
the earth resources survey program. In particular, the 
Range Team was assigned to attempt vegetation type-mapping 
of range and related grazing land. A vegetation type is 
construed as a collection of vegetation of similar, although 
not necessarily identical, species dcxninated by a single 
characteristic species or a snail number of co-dominant 
species. 

To evaluate the utility of ERTS-1 multispectral scanner 
data, the Range Team attempted two levels of classification, 
first-order mapping and second-order mapping. First-order 
mapping was the separation of natural vegetation into the two 
broad categories of woodland and nonwoodland. Second-order 
mapping was the classification of the first-order categories 
into specific vegetation types. The parameters used for 
quantitative analyses were the accuracies of type identifi- 
cation and of areal measurement. The values of these two 
parameters were obtained by analyzing at least two sets of 
data (two dates of the same scene) of two study areas 
selected in the Houston Area Test Site of Texas. The pur- 
pose of performiiig multidate analysis was an attempt to make 
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use of the periodic coverage by ERTS-1 and to investigate 
the possible advantages in data classification and analysis 
gained by a knowledge of the changes %diich occurred between 
different phenological stages. 

Two study areas in the Houston Area Test Site were 
selected for investigation, the Snook site and the San Ber- 
nard site. The Snook site lies astride the boundary between 
the Texas blackland prairie region and the Texas claypan 
savamnah. The features of interest in this site were native 
stands of post oak, bottomland hardwood, n^squite, native 
grassland, and planted pastures of bermuda grasses. The 
San Bernard site is typical of the Gulf Coast marsh region 
along the Texas and Louisiana coasts. Two distinct vegeta- 
tion zones existed in this site, the wet lowland zone of 
marshhay cordgrass and the drier upland zone of gulf cord- 
grass. Both sites offered the opportunity of a first-order 
broad classification, while the San Bernard site provided 
an expected possible second-order classification of vegeta- 
tion types. 

Two classes of techniques were exercised to process the 
multispectral scanner, data which was recorded on film and 
digital tape. They were conventional image interpretation 
techniques and computer-aided data processing techniques. 

In the first class of techniques, multiband ERTS-1 film 
imagery was first displayed on optical and electro-optical 
color film viewers. Color enhancements were created in which 
the features of interest were accented from the background 
scene. Photointerpreters were asked to delineate these 
features (areas that appeared similar) , which were then 
pl€mimetered. Ground-truth information derived from 
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low-altitude aircraft photography was subsequently used to 
check the areal measurements obtained. Regression analyses 
were also made to investigate the significance of the results. 

In the second class of techniques, digital tapes were 
first screened and edited. Processing then took place in 
two different modes, supervised IJVRSYS-type classification, 
and nonsupervised classification using the clustering pro- 
gram ISOCLS. The supervised classification technique 
requires information to train the computer, while the non- 
supervised classification technique requires no training 
information to produce a classification-like clustering 
map, in which the clusters are identified in a postprocessing 
analysis. Both these modes resulted in classification maps 
and measures of classification performance which were eval- 
uated against ground truth. 

Favorable conclusions were reached by exeunining the 
limited amount of quantitative results obtained from both 
types of analyses. 

1.1 CONVENTIONAL DATA PROCESSING 

a. The exsunination of the color enhancements created 

on the film viewers indicated that first-order classification 
of both sites appeared satisfactory. This meant the site 
was separated into wooded rangeland, nonwooded open range- 
land, cropland, water, aud urban area. 

b. The separation of the wet lowland zone of marshhay 
cordgrass from the drier upland zone of gulf cordgrass in 
the San Bernard site was achieved, which was second-order 
feature classification. 



1-4 


c. Regression analyses were conducted to correlate 
the areas interpreted from color enhancements against 
ground truth, A linear model was developed to predict the 
true areal measurement from a measurement obtained from a 
color enhancement created on the available film viewers. 
However, the standard errors of the areal measurements were 
as high as 250 hectares (530 acres) when measuring areas as 
large as 10,500 hectares (22,500 acres). The inordinately 
large error figures were believed to be due to the fuzziness 
in boundaries between the extremely complex vegetation 
zones, as well as to the less-than-ideal quality of the color 
enhcuicements created in this investigation, 

d. Some features, such as cropland, were enhanced 
and accented from the background more vividly in the tem- 
poral enhancements than in signal-date enhancements. 

1.2 COMPUTER DATA PROCESSING 


a. The classification accuracies of training fields 
were high. 3y aggregating second-order fields into first- 
order categories, the first-order classification accuracies 
appeared eve.i higher, with values as high as the mid-ninety 
percentile. 

b. The temporal analyses provided even better classi- 
fication accuracies for the training data as well as less 
confusion between training classes. 

c. The classification accuracies of test fields were 
not as high as those of training fields. Again, first-order 
classification accuracies appeared higher than second-order 
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accuracies. The choice of nonrepresentative test fields was 
believed to be due to the lack of intensive ground truth 
and spot-checked fields. 

d. First-order classification was generally satir- 
factory for both sites; i.e., separation into wooded range- 
land, nonwooded rangeland, cropland, and water. In the San 
Bernard site, the soil moisture content in the marshhay 
cordgrass and that in the gulf cordgrass was sufficiently 
different to permit separation between the two classes. 

e. The clustering results were very useful because 
no information was required to identify spectrally homo- 
geneous areas in the data. Further still, this nonsuper- 
vised classification technique allowed the grouping of 
unknown and/or even siibtle features into spectrally unique 
groups which further analysis could perhaps assign signifi- 
cance. A Scunpling technique to reduce the computation time 
of the clustering process was attempted which resulted in 
positive conclusions, 

f. The clustering results appeared to be very satis- 
factory for first-order classification, and to a certain 
extent for second-order classification. A quantitative 
analysis of the classification accuracies was not made due 
to the lack of intensive ground truth. 

From the limited experience gained by this investiga- 
tion, both the conventional photointerpretation processing 
methods and the computer processing methods had their merits, 
even though the latter class of techniques seemed to be more 
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amenable to systematic quantitative analyses. The following 
positive findings were determined: 

1 • The Range Team demonstrated the utility of the ERTS-1 
MSS data for mapping vegetation types satisfactorily 
into first-order broad categories. Rangeland com- 
prised of woodland and nonwoodland was distinguished 
from nonrangeland comprised of water, urban area, 
and cropland. 

2. Second-order finer classification was also possible, 
depending on the existence of certain not-yet-fully- 
understood characteristics of the features, such as 
soil moisture content. The demonstration of the 
ability to separate the wet lowland zone from the 
drier upland zone, e.g. in the San Bernard site, 
should be of significant value to users interested 
in mapping wetland features in coastal areas. 

3. Temporal analysis, particularly using computer-aided 
techniques in processing digital data, enhanced the 
discrimination between nearly all the vegetation 
types in this investigation. This conclusion 
supported the prediction from mathematical theories. 
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2.0 INTRODUCTION 

2.1 OBJECTIVES 

The ERTS-1 Range Analysis Team conducted a 1-year 
investigation to determine the utility of ERTS-1 remote- 
sens*. r data for vegetation-type mapping of range and 
related grazing land. 

A vegetation type is a collection of vegetation of 
similar, although not necessarily identical, species domi- 
nated by a single characteristic species or a small number 
of co-dominant species. The utility of the ERTS-1 data was 
evaluated using the parameters of the accuracies of type 
identification and of areal measurement. 

A schematic drawing of this investigation and an illus- 
tration of the types of data and products are presented in 
figures 2-1 and 2-2. 


2.2 CLASSIFICATION HIERARCHY 

Rangeland was equated to woodland (stands of trees) 
and nonwoodland (open grassland and brushland) . Many tree 
types and numerous vegetation species could be associated 
with these broad categories. A detailed definition is 
presented in appendix A. 

Two natural levels of classification were demanded in 
this investigation, first-order mapping and second-order 
mapping (table 2-1) . First-order mapping was the categori- 
zation of natural vegetation into the broad categories of 



2-2 







xs 

c m <S) (IS 

o *■> u 

■H .'O -P iji ft) 


4J <a -H 3 x: 

« 4J 6 O CL 



«s 



o a 

Ck 

(U IM 

« o 


xs 

c 


flj 



o ® 


Figure 2--1.— Schematic drawing of the ERTB-1 rangeland investigation. 
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Types of raw data and final products in the rangeland investigatinn. 
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woodland and nonwoodland, as opposed to nonrangeland, which 
was considered to be composed of cropland, water, and urban 
area. The first-order categories were further classified 
into specific classes in the process of second-order mapping. 


TABLE 2-1.- FIRST-ORDER AND SECOND-ORDER 
CLASSIFICATION HIERARCHY^ 


First Order 

Second Order 


Post-oak stands 

/ Woodland 

Mesquite stands 
Bottomland hardwood 


/Open pastures 
1 Abandoned cropland 
■Gulf cordgrass 

1 Nonwoodland 

J Smut cordgrass 
\Marshhay cordgrass 
1 Burned area, previously cord- 


1 grass 

Icultivat-id bermuda fields 
\Wet bottomland 

/ W?ter 

/Deep-sea water 
? Coastal water 

I 

(Inland water, etc. 

Non- / 


rangeland 1 Cropland 

Rice fields, etc. 

lurban area 

Residential, etc. 


^Although not necessarily typical of range features, 
the tabulated second-order features under the class of range- 
land were chosen because they were dominant in the two study 
sites. Attempt was not made in this investigation to classify 
second-order features of nonrangeland, thus they are not 
defined in detail. 
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2.3 SCOPE OF INVESTIGATION 

The investigation area was the Houston Area Test Site 
(HATS) in Texas. Two sites within HATS were selected for 
study. While the results provided insight into the general 
utility of ERTS-type data, they might not be directly appli- 
cable to ecosystems of totally different structures in other 
rangeland areas of the United States. 

Two to three sets of data (two to three dates of the 
same scene) of each of the two study sites were investi- 
gated. The data sets spanned a period of 3-1/2 months, 
from late August to mid-December 1972. This offered the 
opportunity of temporal analyses, which were conducted to 
determine the utility of rhe cyclic coverage of the ERTS-1. 

The term temporal analysis is used in this investiga- 
tion to connotate the analysis of composed data sets. Such 
a composed data set is a unique set of data that is derived 
from considering a feature to be described by a temporal 
signatuj - The term temporal signature is used to construe 
the cha- • .ceristic of a feature across time, with a time- 
span yet to be specified. A simplified example is that of 
deciduous trees: the tree has a 1-year temporal character- 

istic of being light green in spring, richly deep green in 
summer, yellow in autumn, and brown in winter. Temporal 
analyses were conducted to identify the possible advantages 
in data classification and analysis gained by a knowledge 
of the changes which occurred between different phenological 
stages. 
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Since limited manpower and resources prohibited com- 
plete and detailed on-the-spot ground-truth acquisition, 
low-altitude aerial photography was used as a substitute. 

On these photographs, homogeneous areas were delineated, 
from which selected targets were point-checked on the ground. 
The identification of the existing dominant vegetation type 
or types in the targets subsequently produced controlled 
ground-truth information on the study sites. 

The ERTS-1 multispectral scanner (MSS) data in the 
forms of film imagery and digital tape data were used in 
this investigation. These data were acquired from the 
Goddard Space Flight Center (GSFC) . 

The film data were analyzed by conventional photointer- 
pretation techniques and digital tape data were processed by 
computer-aided techniques. A schematic drawing illustrating 
the data flow in this investigation using the two classes of 
methods is shown in figure 2-3. 



riMiMtry 
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Figure 2-3. - Data flow of the ERTS-1 range investigation using conventional and 

computer methods. 
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3.0 STUDY AREAS 

The two areas within HATS selected for study, the 
Snook site and the San Bernard site, are described in detail 
below. 


3.1 SNOOK SITE 

The Snook site is northeast of Somerville, Texas, in 
Burleson County (figures 3-1 and 3-2) . The coordinates of 
the approximate center of the area are latitude 30*26' N. 
longitude 96“27*W. 

The total area of the site is approximately 300 square 
kilometers (116 square miles). The Snook site lies astride 
the boundary between the Texas blackland prairie region and 
the Texas claypan savannah. Thus, this site afforded the 
opportunity to excunine the separability of plant communities 
characteristic of each region in the ERTS-1 imagery. 

The total area of the blackland prairie is about 
51,800 square kilometers (20,000 square miles), most of 
which is in central Texas. Over two-thirds of the area is 
cropland, about one-sixth is pasture, and the remainder is 
low woodland. The woodland areas are characterized by such 
noncommercial species as pecan {Carya sp») , oak (Querous sp,), 
hackberry {Celtis sp,) , and elm {Ulmus sp.) , 

The Texas claypan savannah, about 33,700 square kilo- 
meters (13,000 square miles) in size, is about one-half 
woodland. Most of the wooded area is in post oak (Quercus 
stello.ta Wagenh.) with an understory of yaupon [Ilex 
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vomitoria Ait.) and winged elm (Ulmus alata Michx.). Some 
hardwood forest occurs in the wetter bottomlands. The more 
open areas of post oak have a grass understory of brownseed 
paspalum (Paspalum plicatulum Michx.), silver bluestem 
{Andropogon saooharoides Swartx) , and little bluestem 
{Andropogon Sooparius Michx.). About one-half the area had 
been cleared and cultivated in the past, but most of the 
land has returned to pasture of native grasses and forbs. 

The second-order features (vegetation types) to be 
investigated on the Snook site included native stands of 
post oak, bottomland hardwood, abandoned cropland, mesquite 
[Prosopis sp.), and planted pastures of bermuda grass. A 
first-order breakdown of the ground-truth information on 
the site is shown in figure 3-3. 

3.2 SAN BERNARD SITE 

The San Bernard site is west of the San Bernard River 
outlet on the Gulf of Mexico (figure 3-4), Its approximate 
center is located at latitude 28*42 *N. longitude 95® W. 

The San Bernard site, about 370 square kilometers 
(143 square miles) in area, is typical of the Gulf Coast 
marsh, a 20,700 square kilometer (8,000 square mile) region 
along the Texas and Louisiana coasts. Virtually none of 
the area is farmed, but livestock grazing is a significant 
seasonal use of the drier areas. Hunting and fishing are 
important recreational uses of the region, and trapping 
furbearers is of local economic importance. 
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Figure 3-3.- First-order ground truth of the Snook site. 



Figure 3 ~ 4 .~ Iiifrared aircraft photog'raphic mosaic of tiio 

San Bernard site . 
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The coastal marsh has two fairly distinct vegetation 
zones, marshhay cordgrass (Spartina patens Ait. Muhl.) and 
smooth cordgrass {Spartina alterniflora Loisal.). These 
zones occur on the salt marsh site immediately adjacent to 
the Gulf waters from 0 to 2 inches mean sea level (m.s.l.) 
and are usually inundated during high tide. Adjacent to 
the salt marsh is a slightly higher zone (2 inches to 3 feet 
m.s.l.), dominated by gulf cordgrass {Spartina spartinae 
Trin, Merr.). Each of these two vegetation zones was a 
second -order feature ii. the investigation. A fiist-order 
breakdown cf the ground-truth information on the site is 
shown in figure 3-5, 
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4.0 DATA UTILIZATION 


The ERTS-1 MSS data acquired from the GSFC used for 
analysis were in the form of film imagery and digital tapes. 
Aircraft data were used to create controlled ground-truth 
information. 


4 . 1 ERTS-1 DATA 


The ERTS-1 MSS data were collected in the four spectral 
bands 4 through 7 (equivalent to channels 1 through 4 in 
other references) . Band 4 (channel 1) covers the spectral 
band cf 0.5 to 0.6 micrometer (ym) , band 5 (channel 2) 
covers 0.6 to 0.7 ym, band 6 (channel 3) covers 0.7 to 
0.8 ym, and band 7 (channel 4) covers 0.8 to 1.1 ym. These 
four bands cover most of the visible spectrum and part of 
the reflective infrared spectrum. 

The ERTS-1 data in film and digital tapes which were 
received and processed are tabulated in table 4-1. 

4.2 AIRCRAFT DATA 


Aircraft color and color-infrared photography was used 
as grounr'.-truth information on the distribution cf the vege- 
tation types being studied. The photography was obtained 
from three NASA missions in three photographic scales. 

Refer to table 4-II for the mission numbers, dates on which 
missions took place, the coverage, the altitudes of flight, 
camera systems, film types, and the scales of the imagery. 



TABLE 4-1.- DISPOSITION OF ERTS-1 DATA RECEIVED 


site 

Scene 
Numbe r 

Date 

Date Re< 
For Ana! 

:eived 

vsis 

Conventional 

Processinq 

Computer 

Processinq 

Film 

Tape 

No . 

of ACVP 
Enh. nce- 
ment 

No . 

of MCFV 
Enhance- 
ment 

No . of 
Cluster 
Maps 

No . of 
Classi- 
fication 
Maps 

Snook 

1038-16303 

8/30/72 

12/26/7.; 

12/9/72 

8 

4 

10 

3 

Snook 

1110-16311 

11/10/72 

1/17/73 

- 

3 

5 

- 

- 

Snook 

1128-16311 

11/28/72 

3/2/73 

2/9/73 

1 

- 

13 

3 

Snook 

1146-16311 

12/26/72 

3/20/73 

- 

3 

3 

- 

- 

San Bernard 

1073-16251 

10/4/72 

1/4/73 

1/19/73 

4 

6 

7 

3 

San Bernard 

] 127-16260 

11/27/72 

4/14/73 

4/20/73 

L-:: 

4 

2 

2 
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All delineation and subsequent planinietry of ground- 
truth features were done on positive prints of the photo- 
graphs. First-order type mapping for the two study sites 
was done on 1:120,000 scale photographs. Second-order 
intensive mapping was done on low-altitude 1:20,000 scale 
photographs . 
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5,0 DATA PROCESSING METHODS AND EQUIPMENT 


Two classes of techniques were exercised to process 
the MSS data. Conventional photointerpretation methods 
were employed to process film data, while computer data 
processing methods were employed to process digital tape 
data. 


5.1 CONVENTIONAL DATA PROCESSING 

Conventional photointerpretation data processing of 
film imagery data basically involved two steps, creating 
color enhancements of the imagery on some multichannel 
viewing equipment and analysis of the color enhancements. 

5.1.1 Viewing Equipment 

Two machines were used to create color enhancements, 

the Itek Additive Color Viewer/Printer (ACVP) shown in 

2 

figure 5-1 and the International Imaging Systems (I S) 

Model 2000 Multichannel Film Viewer (MCFV) shown in 
figure 5-2, The former is an optical multiband viewing 
equipment. The latter is an electro-optical imaging 
equipment. 

The Itek ACVP accepts up to four single-band black-and- 
white transparencies and projects them through different 
colored filters onto a common viewing surface. An 8- by 
10-inch film holder fits over the viewing surface to photo- 
graph the color image constructed. 
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Figure 5-1 
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Fiqurc 5-2.- International liiiacjing Systems (I S) Model 2000 Multichannel Film 
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2 

The I S Itodel 2000 MCFV is a more sophisticated type of 

equi(Mnent. A discussion on the operation and capcd>ility of 

2 

the I S machine is found in appendix B and in reference 1, 

5.1.2 Creation of Color Enhancements 

A false color display of the image scene enhances the 
features of interest from the background scene. The 
rationale of creating good colo'. enhancements is discussed 
in appendix B. The two instruments described above were 
used for this process. 

5.1.3 Analysis of Color Enhancements 

After the color enhancements were made, they were inter- 
preter' by three persons not associated with the Range Teeun 
and unfamiliar with the study sites. The features that 
appeared similar were delineated and planiroetered to generate 
area figures. Ground-truth information was compared with 
the areal measurements obtained. A regression analysis was 
then performed to investigate the acctoracy and statistical 
significance of these results. 

5.2 COMPUTER DATA PROCESSING 


The computer processing facilities which were used to 
process the digital I tape data were: 

a. UNIVAC 1108 computers 

b. IBM 360/75 computers 

c. CDC CYBER 73 computers 




Earth EesQurceiS Irtte 
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d. The E?,rth Resources Interactive Processing System 
(ERIPS) , wnich is accessed through two Hazel tine 
terminals to the IBM 360/75 computers in the Real- 
time Computing Complex (RTCC) of the NASA Lyndon 
B, Johnson Space Center (JSC) (figxire 5-3) 

e. The Purdue Terminal, a remote computer terminal 
connected to the IBM 360 computing facilities and 
software packages at Purdue University of Lafayette, 
Indiana 

f . Two Data Analysis Sections (DAS) , on which Computer 
Compatible Tapes (CCT) of multispectral data are 
screened and edited (figure 5-4) 

The software programs which were used for the classi- 
fication and analysis activities are discussed in paragraphs 
5.2.4 and 5.2.5. 

The following describes the five general data processing 
activities: 

1. Tape reformatting 

2. Screening and editing data tapes 

3. Data preprocessing (data registration and 
composition) 

4. Supervised classification (LARSYS) 

5. Nonsupervised classification (ISOCLS clustering) 

5.2.1 Tape Reformatting 

The lack of a common format among all data generation 
and processing systems at the JSC required that certain 
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TABLE 5-1.- REQUIRED TAPE FORMAT CONVERSIONS 


From 

To 

Purpose 

GSFC 

MSDS 

To allow ERTS data to be screened 
on the DAS. 

GSFC 

LARSYS II 

To allow ERTS data to be processed 
by the programs on the UNIVAC 1108 
and the IBM 360/75 computers. 

MSDS 

LARSYS II 

To allow DAS edited tapes to be 
processed by the programs on the 
UNIVAC 1108 and the IBM 360/75 
computers . 

9 -track 

7 -track 

To interface the UNIVAC 1108 
(mostly 7-track) with the 
IBM 360/75 (mostly 9-track) 
computers . 
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data format conversions be performed. These conversions and 
their justifications are shown in table 5-1. 

5.2.2 Screening and Editing 

The ERTS data tapes in multispectral data system (MSDS) 
format were screened on the PMIS DAS. The selected study 
sites were then edited to obtain the desired scale. Edited 
copies of the image at the same scale and at a 2:1 scale 
were usually created. 


5.2.3 Data Preprocessing 

Image-to-image data registration is a requirement for 
temporal analysis and was available on the ERIPS system. 
After two images of the same study site were registered, an 
ERIPS image composition process combined the two 4-channel 
data tapes into an 8-channel tape. 

Another preprocessing activity that was considered was 
atmospheric correction to remove the effect of haze in the 
atmosphere and to accommodate the effect of seasonal changes 
in the Sun angle. This activity was not carried out because 
of the lack of time and the delayed availability of the 
atmospheric correction software. The removal of atmospheric 
effects was also shown not to affect the single-date classi- 
fication of small areas such as the present study sites 
(ref. 2) . 


5.2.4 LARSYS: Supervised Classification 

The ERIPS interactive system was made available to 
Earth Resources users in February 1973. Its software 
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includes a set of programs for multispectral pattern recog- 
nition called LARSYS, The programs were written by the 
Laboratory for Applications of Remote Sensing (LARS) at 
Purdue University, The supervised classification technique 
(ref. 3) involves the following steps: 

a. Select training fields. 

b. Measure the spectral characteristics of the training 
fields by generating statistics of their spectral 
distributions . 

c. Select the channels which will best discriminate 
between various classes. 

d. Classify the test data using the maximum likelihood 
classifier (ref. 4). 

e. Display the classification results on the line 
printer and on the PMIS DAS where color film copies 
can be generated. 

Refer to appendix C for a detailed description of a 
standard LARSYS-type procedure. 

5,2,5 ISCX:lS Clustering: Nonsupervised Classification 

Another method of classifying multispectral data is 
nonsupervised classification or clustering with program 
ISOCLS (refs. 5 and 6). The program groups data points 
that are spectrally similar to create a classification-type 
cluster map. The clusters are then identified with the 
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help of ground truth. This process was carried out in the 
following steps. 

a. The data were input to the Univac 1108 ISOCLS 
program and clustered. 

b. A cluster map was produced on the line printer and 
on color film on the PMIS DAS. 

c. The clusters were identified by checking selected 
areas against ground truth and color coded to 
produce classification maps. This process is 
discussed in more detail in paragraph 6.2.2. 

Appendix D provides a detailed description of the 
procedure to exercise the clustering program. 
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6.0 RESULTS AND DISCUSSION 


This section consists of three main subsections. The 
first subsection presents the processed output products 
using conventional photointerpretation techniques and the 
analyses of these film products. The second subsection 
presents the computer data-processing results of the 
digital tape data. The output products derived by both 
supervised and nonsupervised classification techniques are 
discussed. A general discussion of these two computer 
methods is also included. The third subsection summarizes 
the work that was done and offers suggestions for improving 
future studies. 

The relevant photographic products of this investiga- 
tion are documented in appendixes H, I, and J. Only four 
representative classification maps are illustrated in this 
section. The classification results and cluster interpre- 
tations are also tabulated i*' appendixes I and J, and a sum- 
mary of their analyses is presented in this section, 

6.1 CONVENTIONAL PROCESSING 

6,1.1 Color Enhancements 

Single-date, single-scene color enhancements were 
created on film viewers using different spectral bands of 
the same single-date film imagery on the different channels 
of the viewer. For example, band 4 of the ERTS-1 imagery 
was placed on the blue channel of the viewer, band 5 on the 
green channel, and band 7 on the red channel. Temporal 
color enhancements; i.e., multidate, single-scene enhance- 
ments, were created on the viewer using scenes of different 
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dates on different channels of the viewer. For example, a 
temporal composite of three dates could be made by placing 
band 4 of the second date on the blue channel, band 6 of 
the first date on the green channel, and band 7 of the 
third date on the red channel. 

Typical color enhancements from the ERTS-1 film imagery 
are presented in figure 6-1 for the Snook site and in 
figure 6-2 for the San Bernard site, 

6,1.2 Interpretation 

A good enhancement which offered high contrast between 
different features was not easy to create due to the com- 
plexity of the second-order range features in the Snook 
site. Figure 6-1 in particular shows the fuzziness of the 
boundaries and the gradual transition of color shades 
between features. However, first-order separation of post- 
oak stands, wooded rangeland, bermuda fields, nonwooded 
rangeland, water, and cropland could be detected. The tem- 
poral composite of figure 6-1 indicates that some features, 
such as cropland, are more vividly enhanced and accentuated 
from the background scene in this temporal composite than 
in the three corresponding single-date, single -scene compo- 
sites (figures H-1, H-2, and H-3) , 

Figure 6-2 shows the separation of the wet lowland 
zone of marshhay cordgrass from the drier upland zone of 
gulf cordgrass in the San Bernard site. These were second- 
order features under the category of nonwooded rangeland. 
Woodland, water, and urban areas were also detectable. 




Figure 6-1.- Temporal composite of ERTS-1 data collected 
over the Snook site. (Single-scene enhancement of film 
imagery of August 30, November 10, and December 16, 1972). 
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6.1.3 Areal Measurements and Regression Analysis 

A total of 41 enhancements were produced from the Snook 
and Bernard imagery. Of these, time limitations 
permitted only '' 1 to be interpreted and planimetered 
(figures H-1 througii H-11) . The planimetry was accomplished 
on a Dell Foster data quantizer. After the interpreted 
areas had been determined and tabulated by acreage measure- 
ment, data date, machine used, and the interpreter who 
created the enhancement, this infoinnation was analyzed in 
several ways. The average values from the three delineations 
by three interpreters are presented in tables 6-1 and 6-II. 

A regression analysis was one method used to evaluate 
the derived data. A linear model was chosen represented 
by the following equation; 

Y = A + BX + E 

where Y = areal measurement from ground-truth maps, 

X = areal measurement from planimetering color 
enhancements of ERTS-1 film imagery, 

A,B = coefficients to be estimated, and 
E = error. 

The coefficients A and B were estimated by inputting 
sets of measurements of X with the corresponding ground- 
truth measurements Y into a linear regression program on 
the Univac 1108 computer. A standard error figure e of 
the closeness-of-fit of the prediction model 

Ak A 

Y = A + BX 
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TABLE 6-11.- SAN BERNARD SITE FIELD AREAL MEASUREMENTS 
PLANIMETERED FROM GROUND-TRUTH MAPS AND FROM COLOR 
ENHANCEMENTS OF ERTS-1 FZUl IMAGERY 


Pi«Id 

Mo* 

Feature 

Typo* 

GrouBd'' 

Truth 

weesureuewt 

I^S HeaaurcMnts 
(hectares) 

ACVF HeasureaeBts 
(hectares) 


(hectares) 

Oct • 

Hov 27 

oet 8 

NOT 27 


QC 

1,160.0 

929.8 

- 

788.8 

1,210.8 

5S 

GC 

SC3.6 

588.3 

506.7 

588.6 

801.9 

S3 

OC 

273.7 

363.3 

- 

298.8 

- 

SC 

QC 

68.8 

69.6 

- 

- 

- 

S3, 73 

TK 

saa.c 

- 

- 

588.0 

- 

18,19 

w 

232.0 

- 

215.3 

- 

251.8 

18A 

■ 

28.3 

- 

19.8 

- 

21.5 

19 

H 

«0.1 

- 

19.8 

- 

- 

S9 

H 

12.1 

- 

32.S 

- 

- 

S9,60 

W 

26.3 

- 

25.6 

- 

32.0 

«S 

M 

39.7 

- 

- 

- 

81.0 

1,2 

w 

27.5 

30.8 

- 

38.1 

- 

IS 

w 

100.8 

113.7 

116.8 

285.8 

91.1 

S 

8 

1«.6 

- 

16.3 

11.7 

13.8 

«« 

8 

301.2 

388.1 

- 

- 

- 

2 

T 

3S.2 

60.1 

- 

- 

28.1 

10 

T 

1«.2 

18.6 

- 

- 

- 

16,17, M 

HR 

2.9«2.0 

- 

- 

- 

2,560.9 

M,«7.16 

MR 

2,918.1 

- 

2,792.3 

- 

- 

1 JC 

RB 

2,957.0 

2,719.7 

- 

3,731.7 

- 


*GC « gulf cordgrass, W « wat«r. “ transition (gulf cordgrass to oarsh- 
hay cordgrass), T > wooded, B • burned gulf cordgrass, HR > aarshhay cordgrass. 
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was obtained. The symbol e is actually the statistical 
standard deviation of E . The prediction model states 
that future measurements obtained fr<Mn planiroetering a 
color enhanc^nent under similar conditions as those in the 
range investigation will permit the true areal measurement 
to be predicted by A + BX within an average error of e , 

Table 6-III presents the coefficients A, B of the regres- 
sion models and the standard errors of the estimators for 
the cases when all features were considered, when only 
wooded rangeland was considered, and when only nonwooded 
rangeland was considered. 


6.1.4 Discussion 

The standard errors of the linear regression iiKjdels 
were as high as 250 hectares (530 acres) when measuring 
areas as large as 10,500 hectares (22,500 acres). 

The numerical figures in tables 6-1 and 6-II indi- 
cate the large fluctuations in areal measurements across the 
enhancements from different dates and from different instru- 
ments. While rangeland features normally do not change 
substantially during a few months, such large fluctuations 
could be attributed to the fuzziness of the boundaries 
between the very complex range features. The inordinately 
large error figures were also believed to be due to the 
less-than-ideal quality of the color enhancements created 
in this investigation. 



6-9 


TABLE 6-III.- RESULTS OF LIHEAR REGRESSION ANALYSES ON 
THE AREAL MEASUREMENTS OBTAINED BY PLANIMETERING 
COLOR ENHANCEMENTS OF ERTS-1 FILM IMAGERY 


Site 

Feature 

A 

(acres) 

B 

(dimensionless) 

Standard error 
e (acres) 


All 

8.0 

1.03 

430.5 

Snook 

Woodland 

39.2 

1.04 

529.6 


Non- 

Woodland 

-0.9 

1.00 

311.3 

San 

Bernard 

All 

62.8 

0.96 

513.0 
























6-10 


Time did not allow further analysis into other sources 
of errors, the effects of the dates of the imagery, the 
instruments, and the various channel combinations. Higher- 
order regression models should also be examined. 

6.2 COMPUTER PROCESSING 


The results from supervised and nonsupervised classifi- 
cation are discussed and an overview of the results of 
these two types of computer-processing methods is presented. 
A typical classification map by supervised classification 
techniques and a typical classification-like clustering map 
by nonsupervised classification techniques are presented 
in figures 6-3 and 6-4. 

6.2.1 Supervised Classification 

Two sets of ERTS-1 MSS Snook data and two sets of San 
Bernard data were investigated: the Snook site data were 

acquired August 30 and November 28, 1972; the San Bernard 
data were acquired October 4 and November 27, 1972. 

The two sets of Snook data were registered to each 
other on the ERIPS using Snook-1 as the reference scene. 

The sets then were composed to produce 8-channel data sets. 
The new data tape, Snook-1/2, comprised the first 4 channels 
of Snook-2. The scune procedure was followed for the two 
San Bernard data sets. These two data tapes were then proc- 
essed on the ERIPS, using 4 channels at a time for single- 
scene, single-date analyses, and all 8 channels for temporal 
analyses. 






Wet lowland zone 
Drier upland <on> 
Wooded rargexand 
Water 

Urba area, bare 


green 

yellow 

brown 

blue 

of f “pink 



im. * * if* . 

l£2i2I_J 

Interpretation 
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6. 2. 1.1 Results of single-scene» single-date 
processing . 


6. 2. 1.1.1 Snook-1 and Snook-2: The same training 

fields and test fields were selected simultaneously for 
Snook-1 and Snook-2 by selecting them in the registered 
composed data set Snook-1/2. The coordinates of these 
fields applied to both dates and are tabulated in appendix F. 
The representative fields were (1) wet bottomland, 

(2) bermuda class A, (3) bermuda class B, (4) bottomland 
hardwoods, (5) mesquite, (6) dense oak stand, (7) sparse 
oak stand, (8) abandoned cropland class A, (9) abandoned 
cropland class B, and (10) water (Lake Somerville) . 

Using the standard LARSYS procedure, classification maps 
of Snook-1 and Snook-2 were produced (appendix I, figures I-l 
and 1-2) . The training and test fields are also outlined 
in the color prints. Classification accuracies of the 
training fields and test fields are shown in tables I-I 
through I-IV in appendix I. 

6. 2. 1.1. 2 San Bernard-1 and San Bernard-2: Training 

fields and test fields were selected similarly for SB-1 and 
SB-2 (appendix G) . These fields were representative of 

(1) deep-sea water, (2) coastal water, (3) inland lakes and 
turbid water, (4) marshhay cordgrass, (5) burned gulf cord- 
grass, (6) gulf cordgrass, (7) bare soil, (8) trees, and 
(9) smutgrass, Classifcation maps and tables of classifi- 
cation accuracy are found in figures 1-3 and 1-4 and 
tables I-V through I-VIII. 
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6. 2. 1.2 Results of temporal processing . The composed 
8-channel data of Snook-1/2 and SB-1/2 were classified using 
the same training fields as in the single-scene, single-date 
processing. The classified imagery is shown in figures 1-5 
and 1-6. The classification results are tabulated in 
tables I-IX through I-XII. 

6. 2. 1.3 First-order classification . The first-order 
groupings of the 10 classes in the Snook data which were 
constructed resulted in four groups of features. The 
classes of bermuda A, bermuda B, abandoned cropland A, and 
abandoned cropland B were grouped into the category of 
open range; i.e., nonwooded rangeland. The classes of bot- 
tomland hardwood, mesquite, dense oak stands, and sparse oak 
stands were grouped into the category woodland. Bottomland 
and water remained separate categories. 

These definitions for the Snook data permitted two 
more types of classification accuracy tables to be prepared. 
The first type was the first-order classification accuracy 
table of the training and test fields; i.e., the tabulation 
of fields in the four broad categories. The second type 
was the first-order classification accuracy table of training 
and test categories; i.e,, the tabulation of training and 
test categories in the four broad categories. The Snook-1, 
Snook-2, and Snook-1/2 results are tabulated in tables I-XIII 
through I -XXV. 

6. 2. 1.4 Summary of classification results . Bar charts 
were prepared to further summarize these classification 
accuracies (figures 6-5 through 6-8) , The second-order 



6-15 



Xocj(\o;>c ^ur.ojod 


Figure 6-!3.— Snook data second-order classification accuracy 
bar chart for training fields. 
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Fiold 6-6.— Snook data second-order classification accuracy 

bar chart for test fields. 
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Figure 6-7.— San Bernard data second-order classification accuracy 

bar chart for training fields. 
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Figure 6-8.— San Bernard data second-order classification accuracy 

bar chart for test fields. 
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classification accuracies of the fields and classes for the 
single-scene, single-date data sets were plotted against 
the figures frcan the temporal studies. The overall classi- 
fication accuracies (average classification accuracies) are 
plotted in figures 6-9 and 6-10. 

6. 2. 1.4.1 Snook-1 training field classification: As 

shown in table I-I, the second-order classification accura- 
cies for the training fields were fairly high. From 

tables I-XIII and I-XIV, first-order classification accuracies 
appeared even higher, with values as high as the mid-ninety 
percentile. The lower classification accuracies, such as 
those for mesquite (M) and sparse oak (P) , may be explained 
by the natural similarity between these two classes and 
dense trees with grass under stories. Table I-I indicates 
the classification of a number of data points in M and P 
into dense oak and bermuda fields, all were nonwooded types. 

6. 2. 1.4. 2 Snook-1 test field classification: Test 

field classification accuracies were generally poor 
(table I-II) . This is believed due to the lack of adequate 
ground-truth information which resulted in training and 
test fields being selected which were not the same in all 
respects. In addition, various degrees of tree and grass 
mixtures in some of the categories accounted for low classi- 
fication accuracies. 

6. 2. 1.4. 3 Snook-2 analysis: The same conclusions for 

the training field and test field classification accuracies 
were arrived at concerning the Snook- 2 data. Note that the 
seune training field ^cr dense oak as for Snook-1 has 19 data 
points classified into bottomland hardware (table I-III) , 



percent accuracy percent accuracy 
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Figure 6-9.- Snook data overall classification accuracy. 
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while none were classified into hardware for the earlier 
Snook-1 analysis. Similar confusion occurred between oak 
and hardwood in the classification figures for the category 
hardwood. Both results were probably caused by seasonal 
phenological changes that occurred in the vegetation between 
August and November, but which could not be described 
because of insufficient, continuing ground- truth data. 

6. 2. 1.4. 4 San Bernard analysis: The features of 

interest in the San Bernard site were quite well distin- 
guished from one another. Classification accuracies were 
generally high (shown in tables 1-V through I-VIII) . In the 
SB-1 data, smut cordgrass and gulf cordgrass were confused 
with each other; less confusion occurred in the SB-2 data. 

This difference could again be explained by seasonal changes 
that occurred between the two dates. 

6. 2. 1.5 Temporal analysis of training fields . The 
classification accuracies of the training fields of the two 
study sites were higher than the single-date, single-scene 
figures. This conclusion supports mathematical theory, 
which says that more discrimination exists between two 
statistical classes if more channels of independent infor- 
mation are used to describe the two classes. This trend is 
shown by the abundance of off-diagonal zero entries in 
classifications shown in tables I-IX through I-XII compared 
with tables I-I through I-VIII. The increase in classification 
accuracies was due to the increase in correctly classified 
data points which had been previously confused with other 
classes. 
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6. 2. 1.6 Temporal analysis of test fields . The addi- 
tion of extra channels of independent information has shown 
that in the case of good training fields (where classifica- 
tion accuracies were higher than approximately 80 percent) , 
discrimination between features was higher in temporal 
analyses than in single-date analyses. In other words, a 
data point which belonged to class 1 and which was confused 
as belonging to class 2 in a single-date analysis, even 
though it was not unrepresentative of class 1, would prob- 
ably be classified correctly as class 1 in a temporal 
analysis. On the other hand, a data point of class 1 which 
was not very characteristic of class 1 but which was still 
classified in class 1 during a single-date analysis would 
likely be thresholded out (not classified in) class 1 during 
a temporal analysis. This induction was verified by the 
temporal analysis results of the available test fields. 

Good test fields (classification accuracies above approxi- 
mately 80 percent) had higher classification accuracies, 
while bad test fields had lower classification accuracies in 
temporal studies compared with single-date studies, 

6.2.2 Nonsupervised Classification 

Nonsupervised classification is the grouping of multi- 
spectral data into spectrally homogeneous groups without 
computer training. A few parameters in ISOCLS can be input 
by the user which affect the mode of operation. The three 
parcuneters of particular interest are STDMAX, DLMIN, and 
MAXCLS. 

In program ISOCLS, a class (cluster) is considered to 
be a union of subclasses (subclusters) of nominal sizes. 
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The size of a nominal-size cluster is governed by STDMAX 
(maximum stamdard deviation) . A cluster of a standard 
deviation larger than STDMAX in any one of the channels 
would be divided into two subclusters. And if two clusters 
are spectrally similar and have a measure of similarity 
smaller than DLMIN (minimum distance) , the two clusters are 
combined. The MAXCLS (maximum classes) parameter determines 
the maximum number of clusters generated by the progreun. 

The four different approaches discussed in the following 
sections were exercised on Snook-1 data and the results are 
presented in section 6. 2. 2. 2. Nonsupervised classification 
results on Snook-2, San Bernard-1, and San Bernard-2 are 
presented in section 6. 2. 2. 3. Since no CCT compatible to 
the UNIVAC 1108 computer was available on the temporal 
Snook-1/2 and San Bernard-1/2 data and since the clustering 
progreun on ERIPS was not operational at the time of this 
investigation, nonsupervised temporal analysis was not 
performed. 

6. 2. 2.1 Four approaches to classification . In the 
first method. Program ISOCLS generated the ntwessary nominal- 
sized subclusters, and similar subclusters were combined 
automatically by the DLMIN criterion. This was achieved by 
setting MAXCLS to a large value, settino DLMIN to 3.2 (a 
value suggested in reference 4) in regard to the probability 
of misclassification, and setting STDMAX to values between 
2.5 to 3.5 (for the ERTS-1 data which require discrete values 
from 0 to 127) . This procedure produced an approximate 
second-order classification. 

In the second method. Program ISOCLS grouped the data 
into a predetermined number of spectrally similar groups 
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(MAXCLS) . The DLMIN was set to zero and STDMX was set to 
4.5, which produced a coarser grouping than the previous 
method. Depending upon the potential number of classes and 
subclasses present in the data set, limiting MAXCLS to a 
small value yielded a coarse second-order or first-order 
classification. 

In the third method, the results obtained in the second 
method were manipulated so that spectrally similar groups 
were grouped together. Clusters separated by a distance 
less than 3.2 were grouped together either automatically in 
the "chaining process" of the program, or manually during 
display of the output result. This produced an approximate 
first-order classification based on spectral similarity. 

In the fourth method, the results of ISOCLS clustering 
obtained in either method 1 or 2 were initially interpreted. 
Then, clusters that corresponded to classes of the same or 
similar vegetation types were grouped together to produce a 
first-order classification. 

6. 2. 2. 2 Snook-1 analysis . Five classification maps 
of the Snook-1 data are discussed. 

Figure J-1 in appendix J is the result obtained by 
clustering Snook-1 data according to the first approach 
discussed in the preceding section. Only 10 cf the 25 
clusters obtained were considered major because of the 
noise and clouds present in the data. The interpretatirn 
of these clusters is found in table J-I. 
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Figures J-2 and J-3 are the results of using the second 
approach discussed in section 6. 2. 2.1. The computer run 
corresponding to figure J-2 had MAXCLS set equal to 5, and 
the run corresponding to figure J-3 had MAXCLS equal to 8. 
The interpretation of these clusters is found in tables J-II 
and J-III. 

The results of the computer run corresponding to 
figure J-2 were manipulated as discussed in approach 3 of 
section 6. 2. 2.1 to yield figure J-4, a three-class classi- 
fication map. The three classes were interpreted as vigor- 
ously growing cropland (red) , water (blue) , and other, which 
included rangeland of all vegetation types (yellow) . 

The results of the computer run corresponding to 
figure J-3 were manipulated as discussed in approach 4 of 
section 6. 2. 2.1. This gave the first-order classification 
of figure J-5, where five classes are distinguished as 
vigorously growing cropland (red ) , water (blue) , clouds 
(white), wooded rangeland (brown), and other, ncnwooded 
rangeland (yellow) . 

6. 2. 2. 3 Other analyses . The same approaches were 
applied to obtain second-order classifications of Snook-2, 
San Bernard-1, and San Bernard-2. The classification maps 
and cluster interpretations are shown in figures J-6, J-7, 
j-8, and in tables J-IV, J-V, and J-VI. 

6. 2. 2. 4 Discussion on clustering results . From 
theoretical considerations and past experience with ISOCLS, 
different clustering results would normally be expected 
when different parameter values were input to the program. 
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In all cases, the user should interpret the cluster results 
considering the parameter values used. With regard to the 
possible ways to exercise the ISOCLS program and the corre- 
sponding interpretation of the results, four approaches were 
discussed in sections 6. 2. 2.1 and 6. 2. 2. 2. Nonetheless, 
the identification of the clusters remains as a tedious and 
difficult task. 

Up>on ex 2 «nining the cluster results of the Snook-1 and 
Snook-2 data, the Snook-2 data generally had much lower 
reflecting radiance levels than Snook-1 data. In fact, the 
radiance levels were almost halved in the Snook-2 data 
compared to the Snook-1 data. This difference was believed 
mainly due to the seasonal changes in the Sun angle and the 
advancement in the life cycle of the various vegetation in 
the scene. Whatever their causes, these level differences 
made comparing the two results difficult, since no attempt 
was made to accommodate these effects. The same phenomenon 
occurred in the SB-1 and SB-2 data. 

Because the clustering process is rather time-consuming, 
a sampling approach was deemed desirable if the clustering 
results would not be severely affected. The savings in com- 
putation using a sampling scheme were obvious because in a 
sampling scheme where lines and picture elements (pixels) 
are skipped, essentially a smaller amount of data were 
analyzed. To test this, the computer run corresponding to 
figure J-3 was actually repeated with every third, every 
fifth, and every tenth line and pixel under the same running 
conditions as in the case of no sampling. The cluster centers 
of the four runs (tabulated in table 6-IV) agree closely. 

This suggests that the sampling technique is indeed a viable 
approach to reduce computation time. 
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TABLE 6-IV.- COMPARISON OF CLUSTER CENTERS OBTAINED FROM 
CLUSTERING RUNS WITH VARIOUS SAMPLING SCHEMES 


cluster 

Interpreta- 

tion 

Clouds 

Woods 

Wet 

Bottomland, 

Wetland 

Water 

Grassland, 
Sparse , 
Dry 

Healthy 

Vigorous 

Agricultural 

Fields 

Open Range, 
Grassland , 
Bermuda 

Not So 
Vigorous 
Agricultural 
Fields, 
Grassland 

Sampling 

Scheme 

Scheme 1: no sampliru^ 



Scheme 2i every third 

line and pixel 


Scheme 3: every fifth line and pixel 


Scheme 4 

; every tenth 

line and pixel 


66.3 

29.9 

34.4 

28.0 

38.3 

32.0 

33.2 

32.8 

Channel 1 

64.1 

29.7 

34.8 

29.5 

39.0 

32.5 

31.5 

33.3 


59.9 

29.7 

34.8 

27.6 

38.5 

32.0 

31.6 

32.9 


54.3 

30.5 

35.6 

27.7 

34.7 

32.0 

32.3 

33.7 


65.7 

21.7 

29.8 

17.9 

35.4 

22.3 

26.6 

24.6 

Channel 2 

61.6 

21.7 

29.7 

20.8 

36.8 

22.5 

24.4 

25.0 


56.7 

21,6 

29.5 

17.5 

35.1 

21.9 

24.0 

24.7 


54.4 

22.5 

31.1 

17.1 

29.1 

22.1 

24.5 

26.0 


70.4 

37.0 

33.8 

11.8 

47.4 

72.5 

43.0 

52.6 

Channel 3 

69.3 

36.4 

42.3 

13.6 

47.8 

72.8 

41.9 

51.5 


64.9 

35.3 

41.1 

11.4 

48.4 

71.8 

40.9 

51.2 


58.6 

37.9 

14.5 

8.6 

44.1 

71.0 

25.1 

52.1 


33.6 

20.5 

18.2 

5.5 

24.5 

40.5 

23.0 

29.3 

Channel 4 

32.5 

19.8 

22.1 

10.7 

24.3 

38.4 

22.4 

28.7 


31.2 

19.7 

21.8 

4.4 

25.0 

41.5 

22.2 

28.4 


29.6 

.^1.1 

19.2 

4.1 

23.5 

40.5 

17.0 

28.5 
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6.2.3 General Discussion 

6. 2. 3.1 Supervised classification . The present 
LARSYS-type classification technique is the best technique 
available with the existing state-of-the-art, but requires 
the assumptions that (a) the training fields have Gaussian 
statistical distributions, and (b) the user desires to best 
classify the unknown data into one of the available training 
categories; i.e., to identify the spectral similarity or 
dissimilarity between the unknown data and the known cate- 
gories. If a training field consists of only bermuda grass 
at a praticular growth stage under certain soil and climatic 
conditions, a LARSYS classification only identifies bermuda 
fields and non-bermuda fields in the unknown data which are 
spectrally similar to the bermuda training field. Obviously, 
bermuda fields under other conditions and other fields with 
the same spectral characteristics as the bermuda training 
field would be classified as bermuda. This result is not 
the fault of the procedure. Rather, it is due to the lack 
of discrete spectral distinction within as well as between 
various vegetation types. If the user wishes to classify 
all the bermuda fields which exist in the unknown data, he 
needs to know all the subclasses of bermuda fields in the 
data. Then he must program the computer to recognize all 
these subclasses, a very difficult if not impojcible task. 
Furthermore, the user must devise a way to recognize the 
fields that are classified as bermuda, but which really are 
not bermuda fields. This latter requirement is still without 
a solution, 

6. 2. 3. 2 Nonsupervised classification . To compensate 
for the lack of complete training information, nonsupervised 
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classification is an alternative approach because without 
previous information, the technique groups the unknown data 
into spectrally homogeneous clusters. Although an identifi- 
cation of these clusters would then lead to a classification 
map, the process is tedious and often difficult. Also, as 
discussed in section 6.2,2, the clustering results depend on 
the choice of values of the various input parameters. 

Several ways of judiciously exercising the program and sub- 
sequent interpretation of the results were presented. 

6. 2. 3. 3 Temporal analysis . The bar charts of 
figures 6-5 through 6-10 illustrate the verification of the 
mathematical theory that discrimination between features 
becomes higher in temporal analyses than in single-date 
analyses. This was possible by adding extra channels of 
independent information and the description of a feature by 
a temporal signature. The advantage of temporal analyses 
over single-date aniilyses is obvious from the example of 
deciduous trees. Numerous features appear light green in 
spring, many appear richly green in summer, dozens appear 
yellow in autumn, and scores appear brown in winter. 

However, the 1-year behavior of a feature that appears 
light green in spring, richly green in summer, yellow in 
autumn, and brown in winter, suggests that the feature is 

a deciduous tree. 

6. 2. 3. 4 Test fields versus training fields . The 
limited amount of ground truth available made defining 
enough training fields almost impossible. The sama limita- 
tion also prevented test fields being chosen that belonged 
only to the same general category as the training fields. 
Subsequently, the classification accuracies of test fields 
in this investigation were not high. Future studies using 
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satellite data should acquire more intensive ground truth 
so that test fields and training fields can be chosen of 
the same classes and subclasses. 

6. 2. 3. 5 Refinement of training fields . As suggested 
earlier, training classes should be studied for their homo- 
geneity, especially when more than one field of the same 
training class is selected. This is to determine whether 
subclasses should be defined. However, as mentioned pre- 
viously, such definitions require much more intensive ground 
information than was available in this investigation. 

6 . 2 . 3 . 6 Quantitative analysis and comparison . Que . i- 

titative analysis and comparison were not '..ade betv’een tae 
results of different cluster runs or beeween the results of 
the supervised and nonsupervised classification. The main 
reasons were: '1) The JSC capability was not available in 

time to obtain undistorted classification maps and thereby 
to make areal measurement, comparison, and analysis; (2) no 
attempt was made to accommodate the variation in data values 
in the data sets of different dates; (3) time did not pei .’’it 
further in-depth visual comparison between the various 
classification r» -'ts; and, (4) the statistics of the training 
fields defined i. 2 supervised classification processes 

were accidental! ’ st, which prevented a comparison with the 
cluster statist! • of the nonsupervised classifications. 

6.3 SUMMARY AND SUGGESTIONS 
6.3.1 Summary of Work 

Five ERTS-1 MSS film imagery data sets and four ERTS-1 
MSS digital tape data sets were investigated and discussed. 
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Eleven color enhancements were created and interpreted 
using conventional photointerpretation techniques. Six 
classification maps and eight classification-like cluster 
maps were obtained using computer data-processing techniques. 

6 . 3 . 1 . 1 Output products . 

a. Conventional Data Processing 

1. Single-date enhancements of the August 30, 
November 10, and December 16 Snook data. 

2. Multidate enhancements of the temporal Snook 
scene composed of data of the same three 
dates as (1) . 

3. Single-date enhanc^ents of the October 4 and 
November 27 San Bernard data. 

4. Multidate enhancements of the temporal San 
Bernard scene composed of data of the same two 
dates as (3) . 

b. Computer Data Processing 

1. Classification and cluster maps of the August 30 
and November 28 Snook data. 

2. Classification maps of the temporal Snook scene 
composed of data of the same dates as (1) . 

3. Classification and cluster maps of the October 4 
and November 27 San Bernard data. 

4. Classification maps of the temporal San Bernard 
scene composed of aata of the same dates as (3) . 
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6. 3. 1.2 Analysis work on output products . 

a. Conventional Data Processing 

The 11 color enhancements were delineated, interpreted, 
and planimetcred. The delineation and interpretation were 
performed by three person*? unfamiliar with the study site 
and not associated with tae Range Analysis Team. Plani- 
metering of each delineated area v/as performed three times 
and the results were averaged. The areal measurement 
figures obtained were t 2 d>ulated and analyzed through a 
linear regression analysis. The output products and their 
analyses are discussed in sections 6.1.1 and 6.1.2. 

b. Computer Data Processing 

The classification results were tcibulated and manipu- 
lated to obtain the first-order and second-order classifi- 
cation accuracies. The clusters in the classification-like 
cluster maps were interpreted and also tabulated. The 
output products and their analyses are disucssed in sec- 
tions 6. 2. 1.4, 6. 2. 2. 4, and 6.2.3. 

6.3.2 Suggested Improved Methods 

a. Acquire Intensive Ground Truth 

Information on the growth stage, soil moisture, climatic 
condition, and past history, were necessary in the analysis 
and also for defining accurate test areas versus training 
areas. The failure to collect this type of information 
would prevent comprehensive quantitative analyses and 
comparisons. Also, aircraft data should be acquired con- 
currently with ERTS-1 coverage to construe L accurate con- 
trolled ground-truth information. 
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b. Remove Distortion From Computer Products 

Image- to-ground registration of the ERTS-1 processed 
products and registration tc larger-scale imagery used as 
ground truth should be performed. This will remove the 
inherent distortion that exists in the ERTS-1 products and 
would permit areal measurements amenable to comparison with 
ground-truth information. Additional comprehensive quanti- 
tative analyses could then be performed, 

c. Atmospheric Correction and Signature Extension 

The seasonal migration of the Sun causes variations in 
the reflecting radiance levels of even the saune features 
over a period of time. Removing such an effect or even 
transforming to obtain the absolute reflectivities of features 
would permit a meaningful cc7»parison between different com- 
puter classification and clustering results of data sets 
of different dates. Correction of other atmospheric effects 
such as haze would allow the extension of spectral character- 
istics or signatures and could possibly lead to a signature 
bank of various features. 

d. Further Regression Analysis 

Instead of simply using a linear model, higher-order 
regression models should be considered for analysis of the 
planimetx'ic measurements from the delineated color enhance- 
ments. Also, additional methods need to be evaluated to 
provide reliable inferences from the statistical data being 
investigated, as well as for the prediction model obtained. 
Additional analyses are necessary regarding partitioning of 
the variance about regression, and the use of dummy variables 
to asborb effects such as machine (I^S or Itek) , data dates, 
interpreters, and wavelength bands. 
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e. Emphasize Temporal Analysis 

The limited experience gained in this investigation 
indicated the advantages of temporal analyses o r single- 
date analyses, and future studies should attempt more tem- 
poral analyses. The time span of temporal data sets should 
be extended from 2 or 3 months to perhaps even 1 year. Non- 
supervised classification of temporal data sets should also 
be attempted. 
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7.0 CONCLUSIONS 


The ERTS-1 Range Analysis Team studied two areas in 
the Houston Area Test Site in a 1-year investigation to 
determine the utility of ERTS-1 remote sensor data for 
vegetation type-mapping of range and related grazing land. 

At least two sets of data (two dates of the same scene) of 
the study areas were investigated for the purpose of tem- 
poral analyses. Conventional photointerpretation techniques, 
as well as computer data processing techniques, were used 
for processing MSS film and digital tape data. The following 
conclusions were drawn. 

7.1 CONVENTIONAL DATA PROCESSING 

a. Examination of the color enhancements created on 
the film viewers indicated that first-order classification 
of both sites was satisfactory. Each site was successfully 
separated into wooded rangeland, nonwooded open rangelend, 
cropland, water, and urban area (section 6.1.2). 

b. Separation of the wet lowland zone of marshhay 
cordgrass from the drier upland zone of gulf cordgrass in 
the San Bernard site was achieved. These were second-order 
features (section 6.1.2). 

c. Regression analyses were conducted to correlate 
the areas interpreted from color enhancements against 
ground truth. A linear model was developed that would pre- 
dict the true areal measurement from a measurement obtained 
from a color enhancement created on the available film 
viewers. However, the standard errors of the a^'eal measure- 
ments were as high as 250 hectares (530 acres) when areas 
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as large as 10,500 hectares (22,500 acres) were measured. 

The inordinately large error figures were believed to be 
due to the fuzziness in boundaries between the extremely 
complex vegetation rones, and also to the less-than-ideal 
quality of the colcr enhancements created in this investi- 
gation (section 6.1.4). 

d. Some features, such as cropland, were enhanced 
and accented from the background more vividly in the temporal 
enhcuicements than in the single-date enhancements (sec- 
tion 6.1.2) . 


7.2 COMPUTER DATA PROCESSING 


a. Classification accuracies of training fields were 
high. By aggregating the second-order fields into first- 
order categories, the first-order classification accuraci-:-s 
appeared even higher, with values as high as the mid-ninety 
percentile (section 6. 2. 1.4). 

b. Temporal analyses gave even better classi.fication 
accuracies for the training data. Furthermore, confusion 
between trainir. ^ classes was reduced (sections 6. 2. 1.5 and 
6 . 2 . 3 . 3 ) . 

c. Classification accuracies of test fields were not 
as high as those of training fields. Again, the first-order 
classification accuracies appeared higher than second-order 
classification accuracies. The choice of nonrepresentative 
test fields was believed to be due to the lack of intensive 
ground truth and the insufficiency in the number of spot- 
checked fields (section 6. 2. 1.4). 

d. First-order classification was generally satis- 
factory for both sites; i.e., separation into wooded range- 
land, nonwooded rangeland, cropland, and water. In the San 
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Bernard site, the soil .oisture content in the marshhay 
cordgrass and that in the gulf cordgrass were sufficiently 
different to permit separation between the two classes 
(section 6. 2. 1.5). 

e. Clustering results were very useful. This was 
because no training information was required to identify 
spectrally homogeneous areas in the data. Further still, 
this nonsupervised classification technique allowed the 
grouping of unknown and/or even subtle features into spec- 
trally unique groups which further analysis could perhaps 
assign significance. An attempt to use a sampling tech- 
nique to reduce the computation time of the clustering 
process was also made, resulting in positive conclusions 
(section 6. 2. 2. 4). 

f. Clustering results appeared to be very satisfactory 
for first-order classification, and to a certain extent, for 
second-order classification. A quantitative analysis of the 
classification accuracies was not made due to a lack of 
intensive ground truth (section 6. 2. 2. 4). 

7 . 3 OBSERVATIONS 


An attempt to objectively compare the two classes of 
data processing methods was not made. However, from the 
experience gained in this investigation, computer proc- 
essing methods seem to be more amenable to systematic 
quantitative analysis. This is to be expected when the 
delineation and interpretation of color enhancements is a 
tedious and often difficult task. The results of this 
investigation showed that planiraetric results fluctuated 
considerably, the reason for this being attributed to the 
fuzziness of boundaries between the extremely complex 
features. 
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Finally, it is suggested that future investigators 
should entertain the idea of a possible redefinition of the 
terminologies of features or even the objectives of similar 
studies. From the results of this project, second-order 
features, as defined in the project prior to the actual 
investigation, were often confused with one another, except 
in cases where the soil moisture content in the features 
was vastly different. This fact suggests that perhaps 
features should be defined in terms of their soil moisture 
content and other factors, such as bare soil and land, and 
foliage cover or trees because the ERTS-1 satellite senses 
all data as a mixture of features. Perhaps the detection 
of features should be viewed in a different way as the 
detection of these factors in the features. 

Three positive findings are summarized from this 
effort: 

1. The Range Analysis Team demonstrated the utility of 
of ERTS-1 MSS data to map vegetation types satis- 
factorily into first-order broad categories. 
Rangeland comprising woodland and nonwoodland was 
distinguished from nonrangeland comprising water, 
urban area, and cropland. 

2, Second-order finer classification was also possible, 
depending on the existence of certain not yet fully 
understood characteristics in the features, such 

as soil moisture content. The demonstration of the 
ability to separate the wet lowland zone from the 
drier upland zone; e.g,, in the San Bernard site, 
should be of significant value to users interested 
in mapping wetland features in coastal areas. 
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3. Temporal Analysis, particularly using computer-aided 
techniques in processing digital data, enhanced the 
discrimination between nearly all vegetation types 
in this investigation. This conclusion supports 
the prediction from mathematical theories. 

Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, January 24, 1974 
641-14-07-50-72 
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Traditionally, rangeland has been associated with the 
primary use of domestic livestock grazing. Although this 
is only one of many uses of rangeland, for the moment, 
consider it as the only use, and attempt to define range- 
land, based on the production of forage. Although not 
entirely accurate, assimne for the present that cattle eat 
grass only. Then the definition of rangeland must be 
simple grassland. Immediately, however, the problem of 
defining grassland arises, because virtually all the vege- 
tated surface of the earth has grasses on it, and most, if 
not all, is potential grazing land. Therefore, the equating 
of rangeland to plains and prairies is not sufficient. 

Another criterion which could be used to define range- 
land is historical use. And from the standpoint of live- 
stock production, it is probably a fairly good standard to 
use. If past use is chosen as the criterion, however, 
virtually all the area of the 50 United States must neces- 
sarily be included, with the exception of some of the true 
desert regions of the Southwest. There is little doubt 
that any area even poorly suited to domestic livestock 
grazing has not ^een so used. Note that this definition 
includes vast areas of wooded land, brushland, grassland, 
and every combination of all three. Even today, over 
40 percent of the land area of the 50 states is grazed by 
livestock (U.S, Department of Agriculture, 1962) . Now a 
definition of rangeland is possible because it is accurate, 
if not unique. Range was defined by the Society for Range 
Management as "all land producing native forage for (v/ild 
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or domestic) animal consumption. ..." 
range is virtually all-inclusive ecologically, and that 
fact must be recognized if practical remote sensing appli- 
cations for range management are to be developed. 

Although it was assumed in the foregoing discussion 
that grazing was the only use of rangeland, this concept 
must be discarded if the management of rangeland is con- 
sidered in a modern context. No longer can any federal 
agency view grazing as the only use, or even the best use, 
a priori, for any piece of land. Indeed, any grazing on 
public lands is a current point of controversy, if not 
ecologically, at least politically (ref 8) . 

If consideration is restricted to private rangeland, 
as in virtually all the State of Texas, then of course 
grazing might be defined as the only use. The use of private 
land is the responsibility of the land owner, regardless 
of what the "best" public benefit may be, at least for 
the present. Of more importance to the nation is the admin- 
istration of public rangeland, since it can be used in a 
manner responsive to the needs of the country. Besides 
grazing, other uses of public land which must be considered 
include watershed, wildlife, timber and pulpwood, and num- 
erous forms of recreation. 

To sximmarize the foregoing discussion, it is seen 
that [a.) rangeland is not an ecologically unique type of 
land, (b) grazing is not necessarily the only or best use 
of rangeland, and (c) modern resource management practices 
demand many other usages of rangeland, the best utilization 
of which is of keen interest to mankind. 
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APPENDIX B 

CREATING COLOR ENHANCEMENTS ON THE MODEL 2000 MCFV 

1, Rationale — The rationale behind the operation of 

2 

the I S Multichannel Film Viewer (MCFV) is to create 
color enhancements which tend to 'pull out' those features 
of interest from the rest of the scene. This means contrast 
in colors between features of interest and features of no 
interest. Unless it is the case of a pure black and pure 
white duo-shade enhancement of a two-class classification, 
the degree of contrast and thus, the degree of success is 
rather subjective. Furthermore, enhancements are analyzed 
by human interpreters, who, by definition, are individuals 
and are thus subjective. Therefore, considering the extreme 
versatility of the MCFV, it is obvious that no step-by-step 
operational procedure can be set. Rather, some normal 
procedures to approach the problem are presented in the 
following paragraphs, 

2 

2, Preparation of Data — The I S Model 2000 MCFV is 
an electro-optical instrument designed to accept three 
channels of black-and-white positive or negative multiband 
camera film transparencies. The scanner film gates accept 
roll film of up to 5-inch width or cut film chips. Since 
the bulk MSS ERTS-1 data are in 9-1/2 inch or 70-mm formats, 
appropriate cutting and preparation of film chips are some- 
times necessary. 

3, Registration of Data — To fully utilize the capa- 
bility of the MCFV, three channels are normally used simul- 
taneously, The films previously prepared and mounted on the 
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film transports are registered by adjusting the translation, 
rotational and expansion (z-axis) controls. This is 
achieved by first fixing one channel of data as the refer- 
ence channel. Then, e other channels are registered to 
the first channel by mixing the two channels on-and-off 
and by noticing whether the composite image 'moves' or not. 

A stationary image when mixing is taken in-and-out and means 
a registered image. 

4. Operating the Analog Mode — This mode of operation 
is equivalent to the operation of ordinary additive color 
viewers. Effectively, the MCFV projects the imagery in 
each of the three channels in the colors of blue, green 

and ’-ed, the shades of each color depending on the 'brightness' 
of the scene (i.e., the radiance level of the particular 
picture element) , The intensity of each channel can be 
adjusted also. This mode of operation normal?y gives an 
overall nice image without strong enhancement of specific 
features . 

5, Operating the Digital Mode — This allows specific 
colors to be assigned to various intervals of intensity in 
a single channel. Thus, specific features can easily be 
a'icented against the background. The operator should try 
to assign various highly contrasting colors to different 
features which hopefully have different spectral character- 
istics. This mode of operation can, however, give very 
weird pictures because of the discrete nature of operation. 
Mcst often, the whole scene and the relative locations of 
features in the scene are unnecessarily suppressed. 
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6, Operating the Mixed Mode — This mode of operation 
is the combination of the analog mode and the digital mode. 
This mode is highly desirable because it enA>odies the 
advantages of the two modes of operation. Highly enhanced 

eatures in the midst of a very informational false-color 
background can usually be obtained. 

7. Hardcopy Output — A permanent record of the display 
can be obtained by photographing the display screen with a 
70-mm Hasselblad camera and Ektachrcme color film. The 
settings on the analog processor and digital processor can 
be copied by hand and/or through the use of the online 
recorder. 
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APPENDIX C 

STANDARD LARSYS-TYPE PROCEDURE AT JSC 

1. Preparation of Data — The specific ERTS CCT (the 
one tape out of the four that cover the desired ERTS freun^) 
is selected with the aid of topographic maps available at 
the Project Support Office (PSO) . This ERTS tape is in the 
GSFC format and needs to be reformatted on the DAS to MSDS 
format, during which process scaled edits of tlie study area 
can be made. At -he same time, the editc'd tapes can be 
screened on che DAS, using the MOPS program. If the edits 
are satisfactoiry, they can be reformatted on the DAS to 
LARSYJ II format which is acceptable to the UNIVAC 1108 
computers and the ERIPS interactive system. All these tapes 
are in 9~track 800-bpi formats. 

2. ERIPS Processing — The ERIPS system is a real-time 
interactive system. A complete set of operational procedures 
and tape-handling procedures is available in the user's 
manual (ref. 9) . A successful run on the system normally 
requires 3 to 4 hours. The ERIPS system actually accepts 
9-track 800-bpi CCT's in the LARSYS II format as well as in 
the GSFC ERTS-1 format. These tapes should be made available 
to the IBM 360/75 computer operator in the computer room 
about 10 to 15 minutes before the scheduled ERIPS run-time. 
After appropriate sign-on, the user is ready to perform 
standard LARSYS classification on his data. 

a. Training field selection — In the pattern recogni- 
tion software package of the ERIPS, images can be screened 

Multispectral program on passive microwave imaging 
scanner. 
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in various scales. Areas of interest which contain Known 
features can be delineated through the use of a light nen 
and are designated as training fields of the appropriate 
training classes. These fields can have the shape of a 
concave as well as a convex polyhedron, with a maximum of 
10 sides. 

b. Statistics calculation — The statistics processor 
on the ERIPS C 2 ui be called to calculate the statistical 
mean and covariances of the training fields and training 
classes. These characteristics can be displayed in the 
form of histogreuns and also in the form of reports. 

c. Selection of best cheumels — The 'divergence* 
processor on the ERIPS can be called to select the 'best* 
channels to be used in the classification of the unknown 
data. In this range analysis using 4-channel data in 
single-scene single-date studies, this step was deemed 
unnecessary. 

d. Classification — This can be done on the ERIPS 
itself or on the CYBER 73 computers. In this process, the 
statistical information from the training classes are util- 
ized in the maximum-likelihood classifier on both systems 
to identify the unknown data to the closest (spectrally) 
training class. The ERIPS system generates a bla^k-and -white 
microfiche character classification map as a final product. 

To operate on the CYBER 73 computer, a batch system inter- 
face (BSD tape is generated off the ERIPS. This BSI tape 
contains the unknown data, the statistics information and 
channel information derived from the ERIPS processing. 
Operating in this later mode, a DAS compatible tape of the 
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classified image can be generated, which in turn can be 
displayed in various color combinations on the DAS. Color 
film copies can be made on the DAS. At the same time, the 
CYBER computer also generates a line-printer output of the 
classification map. 

3. Registration and Con^sition of Data — This capa- 
bility is also available on the ERIPS. Two images can be 
registered and composed into a single image with a total 
nundaer of chauinels equal to the sum of the number of channels 
in the two individual images. This was how the temporal 
set Snook- 1/2 and SB-1/2 in this range investigation were 
created. 
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APPENDIX D 

EXERCISING THE JSC CLUSTERING PROGRAM 

1. Preparation of Data — cf appendix C, 

2. Making an ISOCLS Run — The ISOCLS prograun can be 
run in batch modes on the Univac 1108 computers and accepts 
only LARSYS II formatted data CCT, Since the Univac 1108 
computers at JSC operate mostly with 7-track tapes, it is 
good practice to reformat {on the Univac 1108 computer) the 
9-track tapes (generated on the DAS) to 7-track tapes. A 
complete program documentation of ISOCLS is foxind in the 
user's manual (ref, 6). Essentially, the user inputs the 
parameter values he so desires when he submits the card- 
deck to the computer center. The main parameters and their 
recommended settings have been discussed in section 6.2.2. 

A run to cluster 100,000 4-chani.el data roughly takes 

2 minutes per iteration. For this amount of data, the 
Range Analysis Teaun normally considered that 10 to 15 
iterations would produce reasoncibly stabilized results. To 
reduce computer run time, sampling schemes could be used, 
as discussed in section 6.2.4. 

3. Interpretation and Display of Clustering Results — 
The ISOCLS program generates a printout on the line-printer 
and a DAS compatible tape. The clusters can be studiec 
individually on the DAS for their identification. Then a 
color-coded classification map can be generated on the DAS 
on films. It is in this process when clusters are inter- 
preted and identified to known classes by checking clustered 
areas against ground-truth information. 



APPENDIX E 


E-1 


I^S SETTINGS USED ON SPECIFIED FIGURES 

1. SETTINGS FOR THE COLOR ENHT^CEMENT OF FIGURE H-1 
2 

(I S digital multiband composite of August 30, 1972, 
ERTS-1 imagery in bands 5, 6, and 7 of Snook site) 

a. Film 

9-1/2 inch (cut) positive transparencies Photc- 
graphic Technology Division (PTD Copy) frame 
number 1038-1603 

b. Analog Signal Processor Settings 

Channels I and II : Internal, transmit, +, sliced 

Channel III ; Internal, exposure, +, sliced 

c. Photomultiplier 

Channels 1, 2, and 3 ; 1,000 volts (approximate) 

Reference Channel (4) ; 650 volts (approximate) 


1 Channel 1 

Channel 2 

Channel 3 

Offset 

6.64 

7.00 

7.64 

Width 

13.04 

14.15 

14.09 


d. (Digital) Color Contrast Values 


Channel I 

Offset; 

06.64 


Width; 

13.04 

Channel II 

Offset; 

07.00 


Width; 

14.15 

Channel III 

Offset; 

07.64 


Width; 

14.09 
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e. Digital Image Processor Settings 


Slice 

5 

D 

2 

3 

4 



7 


9 

10 

11 

12 

13 

14 

3 

Channel 

Red 

Level 

1 

■ 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Green 

Level 

■ 

■ 















(Band 5) 

Blue 

Level 

















Channel 

Red 

Level 

















2 

Green 

Level 


■ 















(Band 6) 

Blue 

Level 


■ 


1 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■ 

1 

Channel 

Red 

Level 


9 


5 

8 

■ 

■ 

■ 

B 

8 

12 

6 

13 

6 

6 

15 

3 

Green 

Level 


8 


1 

8 

11 

10 

10 

11 

12 

1 

14 

3 

8 

11 

1 

(Band 7 ) 

Blue 

Level 


■ 


■ 

8 

1 

10 

■ 

6 

9 

1 

5 

1 

1 

■ 

1 


•Blank = Neutral Setting 
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2. SETTINGS FOR THE COLOR ENHANCEMENT OF FIGURE H-2 

(I^S digital multiband composite of November 10, 1972, 
ERTS-1 imagery in bands 4, 5, and 7 of Snook site) 

a. Film 

9-1/2 inch (cut) positive transparencies (PTD Copy) 
frame number 1110-16311 

b. Analog Signal Processor Settings 

Channels I, II, and III : Internal, transmit, +, sliced 

c. Photomultiplier 

Channels I, II, and III ; 1,000 volts (approximate) 
Reference Channel (4) ; 600 volts (approximate) 



Channel 1 

Channel 2 

Channel 3 

Offset 

Not 

Recorded 

8.83 

4.91 

Width 

Not 

Recorded 

10.6 

13.92 
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d. Digital Image Processor Settings 


Slice 

B 

B 

2 

3 

4 

5 

6 

■ 

■ 

9 

10 

11 

12 

13 

14 

IS 

Channel 

Red 

Level 


1 

1 

■ 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

(Band 4) 

Green 

Level 


B 

1 

■ 

■ 

■ 

1 

■ 

■ 

■ 

■ 

■ 

1 


1 

1 

Blue 

Level 













1 


8 

1 

Channel 

Red 

Level 


9 

9 

8 

1 

■ 

1 

1 

1 

■ 

■ 

1 

8 


1 

1 

2 

(Band 7) 

Green 

Level 


15 

15 

15 

■ 

■ 

■ 

■ 

■ 

1 

■ 

■ 

1 


1 

1 

Blue 

Level 


15 

1 

■ 

11 

10 

15 

15 

15 

15 

15 

15 

1 


1 

1 

Channel 

Red 

Level 


1 

15 

15 

15 

15 

15 

15 

i 

B 

15 

14 

14 

14 

14 

3 

3 

(Bamd 5) 

Green 

Level 


6 

13 

14 

14 

i 

8 

8 

1 

1 

■ 

1 

1 

1 

1 

1 

Blue 

Level 


1 

1 

1 

1 

1 

1 

14 

14 

1 

1 

1 

5 

5 

1 

1 


*Blank = Neutral Setting 
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3. SETTINGS FOR THE COLOR ENHANCEMENT OF FIGURE H-3 
2 

(I S analog and digital multiband composite of 
December 16, 1972, ERTS-1, imagery in bands 4, 5, and 7 
of Snook site) 

a. Film 

9-1/2 inch (cut) positive transparencies (PTD Copy) 
frame number 1146-16311 

b. Analog Signal Processor Settings 

Channel I (Band 4) ; Internal, exposure, negative, 
normal 

Channel II (Band 5) ; Internal, exposure, +, sliced 
Channel III (Band 7) : Internal, transmit, +, sliced 

c. Photomultiplier 

Channels I, II, and III ; 1,000 volts (approximate) 
Reference Channel (1) ; 600 volts (approximate) 

d. (Digital) Color Contrast Values 

Channel III 

Offset: 8.54 

Width; 10.18 
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e. Digital Image Processor Settings 


Slice 

D 

D 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


Red 

■ 






■ 





■ 

■ 

■ 

■ 

■ 

Channel 

Level 







■ 





■ 

1 

■ 

■ 

■ 

1 

Green 







■ 





■ 


■ 

■ 

■ 

(Band 5) 

Level 

Blue 

Level 

1 

15 

15 

14 

13 

13 

1 

14 

12 

12 

12 

1 

1 

1 

1 

1 


Red 

■ 

■ 

■ 




■ 

■ 






■ 

■ 

■ 

Chcmnel 

Level 

■ 

■ 

■ 




■ 

I 







■ 

■ 

2 

Green 

■ 

■ 

1 

12 

12 

12 

B 


12 

12 



9 


1 

|i 

(Bcuid 6) 

Blue 

Level 

1 

1 

1 



9 

1 

1 





8 

8 

1 

1 

Channel 

Red 

Level 




14 

10 

10 


14 

14 

13 



m 

13 

1 

12 

3 

Green 

Level 



10 

14 

14 








1 




(Bcuid 7) 

Blue 

Level 



13 

13 

13 

. ... J 

13 

14 

10 

8 

1 



1 


1 



♦Blank = Neutral Setting 
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4. SETTINGS FOR THE COLOR ENHANCEMENT OF FIGURE H-4 
2 

(I S analog multiband temporal composite of August 30, 
1972, ERTS-1 imagery band 7 (blue), November 10, 
band 7 (green), and December 16, band 7 (red) of Snook 
site) 

a. Film 

9-1/2 inch (cut) positive transparencies — 
frame number 1038-16303 (Goddard Copy) 
frame number 1110-16311 (PTD Copy) 
frame number 1146-16311 (PTD Copy) 

b. Analog Signal Processor Settings 

Channel I (August 30) 

Band 7 imaged blue 
Internal, transmit, +, sliced 

Channel II (November 10) 

Band 7 imaged green 
Internal, transmit, +, sliced 

Channel III (December 16) 

Band 7 imaged red 

Internal, transmit, +, sliced 

c . Photomultiplier 

Channels I, II, and III; 1,000 volts (approximate) 


Reference Channel (4) ; 650 volts (approximate) 
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5. SETTINGS FOR THE COLOR ENHANCEMENT OF FIGURE H-8 
2 

(I S analog multiband composite of October 4, ERTS-1 
imagery in bands 4 (blue) , 5 (green) , and 7 (red) of 
San Bernard site) 

a. Film 

9-1/2 inch (cut) positive transparencies (PTD Copy) , 
frame number 1073-16251 

b. Analog Signal Processor Settings 
Channel I 

Band 4 — Blue image 
Internal, transmit, +, sliced 

Channel II 

Band 5 — Green image 

Internal, exposure, negative, slice 

Channel III 

Band 7 — Red image 
Internal, transmit, +, slice 

c. Photomultiplier 

Channel I, II, and III; 1,000 volts (approximate) 


Reference Channel (4) ; 600 volts (approximate) 
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6. SETTINGS FOR THE COLOR ENHANCEMENT OF FIGURE H-9 
2 

(I S digital multiband composite of November 27 , 1972, 
ERTS**1 imagery in beuids 4, 5, and 7 of San Bernard site) 

a. Film 

9-1/2 inch (cut) positive transparencies (PTD Copy) , 
frame number 1127-16260 

b. No analog or digital processor settings recorded 
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7. SETTINGS FOR THE COLOR ENHANCEMENT OF FIGURE H-10 
2 

(I S digital multiband composite of November 27, 1972, 
ERTS-1 imagery in bands 4, 5, and 7 of San Bernard site) 

a. Film 

9-1/2 inch (cut) positive transparencies (PTD Copy) , 
ERTS frcune number 1127-16260 

b. No analog or digital processor settings recorded 



APPEUDIX I 
lOCAtZOM or SNOOK SI 


1. Coordinates of the training fields of i 
reference to the Snook-1 scene (1038' 


Field 


First 

Coor. 

(line 

sanple) 

Second 

Coor. 

(line 

sample) 

Bottondand 

B 

932 

1984 

941 

1991 

Berauda, 
Class A 

C 

870 

2017 

867 

2024 

Bermuda, 
Class B 

0 

926 

2073 

918 

2090 

BottoBiland 

Hardwood 

H 

891 

1888 

897 

1905 

Mesquite 

M 

774 

2093 

771 

2098 

Oak, dense 

0 

922 

2025 

925 

2033 

Oak, sparse 

P 

1 

i 

836 

2056 

843 

2060 

Abandoned 
Croplc. .d. 
Class A 

H 

825 

2065 

821 

2067 

Abandoned 
Cropland, 
Class B 

X 

827 

2067 

826 

2072 

Water 

z 

991 

2084 

990 

2108 
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TE FIELDS 


:he Snook study site with 
•16303-MB-3) 


Third 

Coor. 

(line 

sample) 

Fourth 

Coor. 

(line 

sample) 

Fifth 

Coor. 

(line 

sample) 

Sixth 

Coor. 

(line 

samiple) 

942 

1979 

949 

1976 

935 

1964 

935 

1976 

862 

2022 

866 

2013 



914 

2088 

918 

2083 

921 

2070 


889 

1910 

884 

1893 



770 

2103 

762 

2097 

766 

2096 

766 

2089 

918 

2038 

904 

2032 

912 

2024 


841 

2065 

836 

2066 

838 

2062 

836 

2061 

821 

2073 

825 

2073 

826 

2067 


829 

2072 

1 

832 

2072 

832 

2071 



976 

2110 


974 

2085 
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2. Coordinates of the test fields of the Snook study site with 
reference to the Snook-1 scene ( 1038-16 303-MB-3) 


■ 

First 

Coor. 

(line 

sample) 

Second 

Coor. 

(line 

sample) 

Third 

Coor- 

(line 

sairple) 

Fourth 

Coor. 

(line 

sample) 

Fifth 

Coor. 

(I'ne 

saiup^e) 

Sixth 

Coor. 

(line 

sample) 

Cl 

869 

2016 

862 

2021 

859 

2026 




C2 

866 

2068 

862 

2074 

855 

2067 

861 

2065 



D 

9*1 

7055 

944 

2062 

942 

2066 

939 

2061 

939 

2057 


H 

898 

1865 

906 

1889 

901 

1890 

894 

1875 



0 

861 

2000 

864 

2006 

861 

2014 

853 

2010 

856 

2006 

857 

2007 

X 

830 

2065 

825 

2061 

825 

2051 

829 

2050 



z 

981 

2101 

981 

2114 

973 

2113 

i 972 
2101 

980 { 

2102 1 



Field 
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APPENDIX ; 

IjOCATION of SAN BERtUUtO SITE FIELDS 


I. Coordinates of the training fields of the San Bernard study site 
with reference to the SB-l scene (1073-162S1-MB-3) 


Field 

■ 

First 

Coor. 

(line 

sample) 

Second 

Coor. 

(line 

sample) 

Third 

Coor. 

(line 

sample) 

Fourth 

Coor. 

(_ine 

sample) 

Fifth 

Coor. 

(line 

sample) 

Sixth 
Coor. 
(line 
sample ) 

Burned 

Gulf Cordgrass 

B 

1126 

1911 

1130 

1920 

1115 

1922 

1115 

1912 



Bare Soil 

C 

1130 

1730 

1128 

1726 

1129 

1724 

1126 

1710 

1127 

1711 

1131 

1724 

Gulf 

Cordgrass 

G 

1101 

1925 

1093 

1920 

1100 

1896 

1107 

1900 



Karshhay 

Cordgrass 

Ml 

1139 

1886 

1133 

1898 

1128 

1898 

1130 

1888 



Marshhay 

Cordgrass 

H2 

1121 

1850 

1122 

1864 

1112 

1867 

1111 

1843 



Oak 

0 

1051 

1922 

1047 

1936 

1040 

1936 

1039 

1919 



Smutgrass 

S 

1001 

IcoS 

986 

1 Ort » 

986 

996 

* 



Coastal 

Water 

X 

1200 

1847 

1205 

1854 

1202 

1862 

1195 

1866 

1182 

1851 


Inland 

Water 

Y 

1153 

1862 

1159 

1870 

1160 

1880 

1152 

1881 

1146 

1868 

1151 

1862 

Deep-sei 

Water 

Z 

1210 

1669 

1210 

1701 

1183 

1701 

1183 

1670 




2. Coordinates of the test fields of the San Bernard study site with reference 
to the SB-1 scene (1073-16251-MB-3) 


Field 

Symbol 

First 

Coor. 

(lire 

sample) 

Second 
Coor . 
(line 
sample) 

Third 

Coor. 

(line 

sample) 

Fourth 
Coor . 
(line 
sample) 

Fifth 
Coor. 
(line 
sam. Ir) 

Sixth 

Coor. 

(line 

sample) 

Burned 

B 



1199 

1205 



Gulf Cordgrass 


1927 

Is9 

1950 

1926 



Bare Soil 

C 

1189 

1179 

1180 

1191 





1874 

1859 

1838 

1874 



Gulf 

G 

1068 

1070 

1056 

1056 



Cordgrass 


1890 

1905 

1905 

1896 



Marshhay 

M 

1138 

1133 

1126 

1130 

1139 


Cordgrass 


1885 

1896 

1896 

1887 

1884 
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APPENDIX li 


COLOR ENHANCEMENTS OF ERTS-1 FILM IMAGERY 


PRODUCED ON THE IS MODEL 2000 «CF¥ 


AND ON THE ITEK ACVP 


2 

Figure H-l,- I S digital laultiband composite of bands 5, 
6, and 7 {ERTS-1 film imagery of Snook site acquired 
August 30, 1972). 

ORIGINAL FAGE E5 
1^’ POOR OyAi-OT 



















Figure H-4.- I^S analog multibami tempoxai composite of 
ERTS-1 film imagery of Snook site - imagery of August 30, 
1972, band 7 (blue); November 10, 1972, band 7 (green); 
December 16, 1972, band 7 (red). 









gure H-6.- ACVP multiband composite of Snook site (ERTS-1 film imagery acquired 
Decorabor 16, 1972 - band 4, blue; band 5, cfrecn; band 7, red). 




H-? 



Fiqure H-7.- ACVP maltiband temporal composite of ERTS-1 film intiigery of Snook 
site (band 7, red, acquired August 30, 1972; band 7, green, acquired 
November 10, 1972) . 








Met lowland zone of marshhay 
cordgrass 


Gulf of Mexico 


Figure H“8.~ IS analog multiband composite of San Bernard 
site (ERTS-1 film imagery acquired October 4 , 1972 ; band 
4, blue; band 5, green; and: band 1 , red). 






Clouds 


Gulf of Mexico 


2 

Figure H~9.- I S digital multiband composite of bands 4 , 5 , 
and 7 (ERTS-1 fxlro imagery of San Bernard site acquired 
November 27 , 1^72) . 






-1 



|Wel: .lowland of marshhay 
cordgrass 


Wet lowland zone of marshhay 
cx?rdgrass 


Gulf of .’’exico 


2 

iigure 1 S digital nultibanr composite oC bands 4, 5, 

and 7 (EETS-1 film imagery of San Bernard site acquired 
Hovember 27, 1972 1 . 






H-ll 


MASA S-73-31456 



Wet lowland zone of marsWiay 
coragrass 


Gulf of Mejeicoi 


H-ll.- ACVP multiband composite of band 4, blue? 

. 5# green; and band 7 , red CERTS- 1 film imagery of 
Bernard site acquired October 4, 1972) . 
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CYBER 73 LARSYS classification of Snook-1 data 
vxith 2-percent thresholding. 
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TABLE I- I I 


SNOOK-1 SECOND-ORDER CLASSIFICATION PERFORMANCE SUMMARY OF TEST FIELDS 


Claaa 

Do. of 
Samplea 

Percent 

Correct 

Diatributlon of Samolea 


■■■■ 


C 

B 

R 

M 


P 


X 

Z 


Cl 

24 

58.3 

BS 

14 

0 

0 

9 

0 

1 

0 

0 

0 

0 

C2 

60 

78.3 


47 

0 

0 

7 

4 

1 

0 

0 

0 

1 

D 

39 

43.6 

2 

3 

17 

1 


1 

2 

6 

7 

0 

0 

H 

143 

72.0 

0 

0 

0 

103 


10 

19 

9 

0 

0 

2 

0 

87 

44.8 

0 

0 

0 

21 

0 

39 

11 

0 

1 

0 

13 

X 

67 

71.6 

0 

1 

0 

0 

0 

0 

0 

0 

48 

0 

18 

7. 

119 

89.9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

107 

12 

Whole 

Teat 

Area 

122,500 


3pS66 

20,420 

6,784 

6,908 

15,B57 

10,926 

14,031 

7,549 

9,492 

1,787 

15,170 


Overall Clastlf ication Accuracy - 56.5 percent 


Cl ■ Bermuda, claae A, field 1 
H - Bottomland hardwood 
Z ■ Water 


C2 - Bermuda, elaae A, field 2 
O ■ Oak 

Three. ■ Threeholded 


D • Bermuda, claaa B 
X ■ Abandoned cropland, claaa B 

































) !* y *< i » ""^ ^^^ ♦ -’ < 

*’■ i^L -^..yti-.y^ 

Ttpipr^ 



Ket bottom! and (turquoise 
Bermuda, class A (yellow) 

Bermuda, class B (light green) 

Bottomland hardwood (light brown) 
Mesquite (orange) 

Dense oak (dark brown) 

Sparse oak (medium brown) 

Abandoned cropland, class A (light gray) 
Abandoned cropland, class B (dark green) 
Water (blue) 

Null class (white) 


Figure 1-2, 


CYBER 73 LARSYS classification of Snook-2 data 
with 2-percent thresholding . 
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TABLE I-III.- SNOOK-2 SECOND-ORDER CLASSIFICATION PERFORMANCE SUMMARY OP TRAINING FIELDS 


Class 

No. of 
Samples 

Percent 

Correct 



IHHB 

B 


0 

H 

ft 

0 


w 

X 

1 

Eiam 

B 

219 

86.0 

189 

■5 

0 

11 

1 

0 

M 

1 

6 

0 

2 

c 

56 

93.0 

0 

52 

4 

0 

0 

0 

K 

0 

0 

0 

0 

D 

94 

67.0 

0 

6 

63 

0 

5 

0 


2 

15 

0 

3 

H 

170 

79.0 

6 

0 


135 

1 

BH 

17 

0 

0 

0 

2 

M 

93 

55.0 

3 

0 


0 

SI 

mm 

14 

4 

15 

0 

0 

0 

250 

90.0 

0 

0 


19 

0 

226 

3 

0 

0 

0 

2 

P 

55 

47.0 

3 

0 

u 

19 

4 

2 

26 

0 

1 

0 

0 

w 

42 

76.0 

0 

0 

3 

C 

1 

0 

0 

32 

5 

0 

1 

X 

25 

64.0 

1 

0 

1 

0 

3 

0 

0 

4 

16 

0 

0 

LJ 

399 

97.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

388 

11 


Overall Clasilfication Accuracy “ 7S.4 percent 

B • Bottomland C • Bermuda > class A D • oarmuda, class B 

H • Bottomland hardwood M ■ Heaquita O - Oak (denaa) 

P • Oak (sparse) M • Abandoned cropland, class A X - Abandoned cropland, class B 

Thraa. ■ Thraaholdad 


Z • Water 






































TABLE I-IV 


SNOOK-2 SECOND-ORDER CLASSIFICATION PERFORMANCE SUMMARY OF TEST FIELDS 


Class 

Ho. of 
Samplei 


Distribution of Samples 


■■■ 

B 

c 

D 

H 

M 

0 



X 

z 

Kirrm« 

Cl 

24 


0 

5 

17 

X 

1 

0 

0 

0 

0 

0 

0 

C2 

60 

46.7 

0 

28 

6 

2 

4 

5 

2 

0 

0 

0 

11 

D 

39 

3«.0 

0 

5 

14 

0 

6 

0 

1 

1 

3 

0 

9 

H 

143 

62.0 

17 

0 

0 

89 

2 

8 

22 

0 


0 

5 

0 

B7 

26.0 

2 

IS 

2 

32 

5 

23 

8 

0 

sss 

0 

15 

X 



67 

9.0 

0 

Si 

12 

0 

8 

0 

3 

27 


0 

11 

2 

119 

98.0 

0 


0 

0 


mmiiQui 

0 

0 

0 

117 

2 

Whole 

Test 

Area 

122,500 


6,579 

5,538 

15,846 

12,922 



11,052 



6,268 

10,025 

2,796 

27,519 


Overall Classification Accuracy • 42.7 

Cl - Bernuda, class A, field 1 C2 - Bermuda, class A, field 2 

H • Bottomland hardvood O • Oaks 

E * Water * ■ Abandoned cropland, class B 

Threa. - Thresholded D - Bermuda, class B 























































Deep-sea water (dark blue} 

Coastal water (turquoise) 

Inland water (light blue) 

Marshhay cordgrass (light green) 
burned gulf cordgrass (dark green) 


(8) Bare soil (off-pink) 

(9) Oak trees (dark brown) 
noi Null class (white) 


3.- CYBER 73 LARSYS classif ication of San Bernard-1 






I 

CO 


TABLE I-V.- SB-1 SECOND-ORDER CLASSIFICATION PERFORMANCE SUMMARY OF TRAINING FIELDS 


Class 

No. of 
Samples 

Percent 

Correct 

Distribution of Samoles 1 

B 

c 

G 

R 


S— T— 

Y 

Z 

Thres. 

B 

148 

95.3 

141 

0 

0 

6 

0 

0 

0 

0 

0 

1 

c 

36 

94.4 

0 

34 

0 

0 

0 

1 

1 

0 

0 

0 

G 

229 

77.7 

0 

1 

178 

0 

0 

49 

0 

0 

0 

1 

Ml 

85 

95.5 

1 

0 

0 

81 

0 

0 

0 

0 

0 

3 

M2 

201 

81.1 

1 

0 

0 

197 

0 

0 

0 

0 

0 

1 

0 

165 

96.4 

0 

1 

0 

0 

159 

0 

0 

0 

0 

5 

S 

230 

77.4 

0 

0 

48 

1 

0 

178 

0 

0 

0 

3 

X 

263 

93.2 

0 

0 

0 

0 

0 

0 

245 

0 

0 

18 

Y 

187 

97.9 

0 

1 

0 

0 

0 

0 

0 

183 

0 

3 

Z 

910 

97.5 

0 

0 

0 

_J 

0 

0 

0 

0 

887 

23 


Overall Classification Accuracy • 92.3 

B « Burned gulf cordgrass C « Bare soil G * Gulf cordgrass 

Ml a Marshhay cordgrass, field 1 M2 " Marshhay cordgrass, field 2 0 • Oak 

S - Smutgrass X a Coastal water Y - Inland water 


z a Deep-sea water 
Thres. a Thresholded 












































i 



Deep-sea water {dark blue) 

Coastal water (turquoise) 

Inland water (light blue) 

Marshhay cordgrass (light green) 
bui-ned gulf cordgrass (dark green) 
Gulf cordgrass (dull yellow) 
Smutgrass (light yellow) 

Bare soil (off-pink) 

Oak trees (dark brown) 

Null class (white) 


CYBER 71 LARSYS classification of San Bernard-2 
data with 2 -percent thresholding. 



PERFORhIANCE SUMMARY OF TRAINING FIELDS 


1-11 
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( 1 ) 
( 2 | 
13) 
HI 
(51 
C6| 
(7) 
{81 
S (9) 


*et bottomland 
Bermuda, class A (yellow) 

Bermuda, class B (light green) 


(4 1 Bottomland har ' 

(5| Mesquite (or an 

(6) Dense oak (dar 

(7) Sparse oak (me 

(8) Abandoned cropland, class A (light 

(9) Abandoned cropland, class B (dark < 

(10) Water (blue) 

(11) Null class (white) 


5.- CYBER 73 LARSYS classification of Snook-1/2 
temporal data with 2-percent thresholding. 
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ib) burned gulf cordgrass (dark green 
(6) Gulf cordgrass (dull yellow) 
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TABLE I-XIII.- SNOOK-1 FIRST-ORDER CLASSIFICATION PERFORMANCE SUMMARY OF TRAINING FIELDS 


Class 

No. of 
Samples 

Percent 

Correct 

Distribution of Samoles 1 

BOTTOMLAND 

OPEN RANGE 

WATER 

WOODLAND 

THRESHOLDED 

B 

219 

90.3 

197 

17 

0 

0 

5 

C 

56 

93.0 

0 

52 

0 

3 

1 

D 

94 

95.0 

3 

89 

0 

2 

0 

H 

170 

95.0 

0 

8 

0 

161 

1 

M 

93 

90.0 

0 

10 

0 

83 

0 

0 

250 

99.0 

0 

0 

0 

247 

3 

P 

55 

85.0 

0 

8 

0 

47 

0 

w 

42 

81.0 

0 

34 

0 

8 

0 

X 

25 

96.0 

0 

24 

0 

1 

0 

z 

399 

95.0 

0 

0 

379 

0 

20 


Overall Classification Accuracy " 91.9 percent 
BOTTOMLAND > Bottomland (B) 

OPEN RANGE " Bermuda, class A (C) ; Bermuda, class B (D) ; Abandoned cropland, class A (H) ; 
Abandoned cropland , class B (X) 

WOODLAND > Mesquite (M) ; Bottomland hardwood (H) ; Dense oak (0) ; Sparse oak (P) 

WATER • Water (Z) 




TABLE I-XIV.- SNOOK-1 FIRST-ORDER CLASSIFICATION PERFORMANCE 


SUMMARY OF TRAINING CATEGORIES 


CATEGORY 

NO. OF 
SAMPLES 

% 

CORRECT 

DISTRIBUTION OF SAMPLES 1 

BOTTOMLAND 

OPEN RANGE 

WATER 

WOODLAND 

THRESHOLDED 

BOTTOMLAND 

219 

90.0 

197 

17 

0 

0 

5 

OPEN RANGE 

217 

' 1 

92,0 

3 

199 

0 

14 

1 

WATER 

399 

95.0 

0 

0 

379 

0 

20 

WOODLAND 

568 

95.0 

0 

26 

0 

538 

4 


Overall Classification Accuracy = 93.0 percent 


BOTTOMLAND : 
OPEN RANGE; 

WOODLAND ; 


Bottomland (B) 

Bermuda, class A (C) ; Bermuda, class B (D) ; Abandoned cropland, 
class A (W) ; Abandoned cropland, class B (X) 

Mesquite (M) ; Bottomland hardwood (H) ; Dense oak (0) ; Sparse 
oak (P) 


WATER 


Water (Z) 




































TABLE I-XV.- SNOOK-1 FIRST-ORDER CLASSIFICATION PERFORMANCE 


O ^ 


►tJ o 

o g 
o c 



SUMMARY OF TEST FIELDS 


CLASS 

NO. OF 
SAMPLES 

% 

CORRECT 

DISTRIBUTION OF SAMPLES 

BOTTOMLAND 


WATER 


THRESHOLDED 

Cl 

24 

58.0 

0 

14 

0 

10 

0 

C2 

60 

78.0 

0 

47 

0 

12 

1 

D 

39 

85.0 

2 

33 

0 

4 

0 

H 

143 

92.0 

0 

9 

0 

132 

2 

0 

87 

84.0 

0 

1 

0 

73 

13 

X 

67 

73.0 

0 

49 

0 

0 

18 

z 

247 

90.0 

0 

14 

0 

218 

IS 

WHOLE 

TEST 

AREA 

122,500 


3,566 

44,245 

1,767 

47,724 

2,578 


Oveir.ll Classification Accuracy « 80.0 percant 


BOTTOMLAND: 
OPEN RANGE: 

WOODLAND : 
WATER : 


Bottomland (B) 

Bermuda, class A (C) > Bermuda, class U (0) t Abandoned cropland, class A (W) ; Abandoned 
cropland, class B (X) 

Mesquite (M) ; Bottomland hardwood (H) ; Dense oak (O) ; Sparse oak (P) 

Water (Z) 















































table 1“Xvi.- snook-1 first-order classification performance 

SUMMARY OF TEST CATEGORIES 


CATEGORY 

SAMPLES 

% 

CORRECT 

DISTRIBUfl 

[ON OF Si 

MIPLES 

BOTTOMLAND 

OPEN RANGE 

WATER 

WOODLAND 

THRESHOLDED 

OPEN RANGE 

245 

71.0 

2 

173 

0 

49 

21 

WATER 

119 

90.0 

0 

0 

107 

0 

12 

WOODLAND 

247 

90.0 

0 

14 

0 

218 

15 

WHOLE 

TEST 

AREA 

122,500 


3,566 

44,245 

1,787 

47,724 

25,178 


Overall Classification Accuracy = 83.7 percent 


BOTTOMLAND : 
OPEN RANGE; 

WOODLAND : 

WATER ; 


Bottomland (B) 

Bermuda, class A (C) ; Bermuda, class B (D) ; Abandoned cropland, 
class A (W) ; Abandoned cropland, class B (X) 

Mesquite (M) ; Bottomland hardwood (H); Dense oak (O) ; Sparse 
oak (P) 

Water (Z) 




TABLE I -XVI I 


SNCOK-2 FIRST-ORDER CLASSIFICATION PERFORMANCE 


CLASS 

NO. OF 
SAMPLES 

% 

CORRECT 

DISTRIBUTION OF SAMPLES 

BOTTOMLAND 

OPEN RANGE 

WATER 

WOODLAND 

THRESHOLDED 

B 

219 

86.0 

189 

7 

0 

21 

2 

C 

5f 

100.0 

0 

56 

0 

0 

0 

D 

94 

91.0 

0 

86 

0 

5 

3 

H 

170 

95.0 

6 

0 

0 

162 

2 

M 

93 

70.0 

3 

25 

0 

65 

0 

■ 

0 

250 

99.0 

0 

0 

0 

248 

2 

B 

55 

93.0 

3 

1 

0 

51 

0 

1 

W 

42 

95.0 

0 

40 

0 

1 

1 


25 

84.0 

1 

21 

0 

3 

0 


399 

97.0 

0 

0 

388 

0 

11 


Overall Classification Accuracy - 91.0 percent 
BOTTOHIAND: Bottomland <B) 

OPEN RANGE: Bermuda, class A (C) ; Bermuda, class B (D) ; Abandoned cropland, class A (W) ; Abandoned 

cropland, class B (X) 

WOODLAND : Mesquite (M) ; Bottomland hardwood (H)t Dense oak (0) i Sparse oak (P) 

Water (Z) 


WATER 


I 

















































TABLE I-XVIII 


SNOOK-2 FIRST-ORDER CLASSIFICATION PERFORMANCE 


SUMMARY OF TRAINING CATEGORIES 



NO, OF 


DISTRIBUTION OF SAMPLES 

CATEGORY 

SAMPLES 

CORRECT 

BOTTOMLAND 

OPEN RANGE 

WATER 

V700DLAND 

THRESHOLDED 

BOTTOMLAND 

219 

86.0 

189 

7 

0 

21 

2 

OPEN RANGE 

217 

94.0 


203 

0 

9 

4 

WATER 

399 

97.0 

0 

0 

388 

0 

11 

WOODLAND 

568 

93.0 

12 

26 

0 

526 

4 


Overall Classification Accuracy = 92.? percent 


BOTTOMLAND: Bottomland (B) 


OPEN RANGE: 


Bermuda, class A (C) ; Bermuda, class B (D) ; Abandoned cropland, 
class A (W) ; Tdaandoned cropland, class B (X) 


WOODLAND 


Mesquite (M) ; Bottomland hardwood (H) ; Dense oak (O) ; Sparse 
oak (P) 


WATER 


Water (Z) 



















































































TABLE I-XX.- SNOOK-2 FIRST-ORDER CLASSIFICATION PERFORMANCE 


SUMMARY OF TEST CATEGORIES 


CATEGORY 

NO. OF 
SAMPLES 

% 

CORRECT 

riSTRIBUTION OF SAMPLES 

BOTTOMLAND 

OPEN RANGE 

WATER 

WOODLAND 

THRESHOLDEO 

OPEN RANGE 

245 

61.0 

0 

149 

0 

51 

45 

WATER 

119 

98.0 

0 

0 

117 

0 

2 

WOODLAND 

247 

81.0 

20 

5 

0 

199 

23 

WHOLE 

TEST 

ARIA 

122,500 


6,579 

37,657 

2,796 

47,949 

27,519 


Overall Classification Accuracy = 80.0 percent 


BOTTOMLAND; Bottomland (B) 


OPCN RANGE: 


Bermuda, class A (C) ; Bermuda, class B (D) ; Abandoned cropland, 
class A (W) ; Abandoned cropland, class B (X) 


WOODLAND 


Mesquite (M) ; Bottomland hardwood (H) ; Dense oak (O) ; Sparse 
oak (P) 


WATER 


Water (Z) 




TABLE I-XXI.- SNOOK-1 /2 FIRST-ORDER CLASSIFICATION PERFORMANCE 
SUMMARY OF TRAINING CATEGORIES 



NO. OF 

% 

DISTRIBUTION OP SAMPLES 

CATEGORY 

SAMPLES 

CORRECT 

BOTTOMLAND 

OPEN RANGE 

WATER 

WOODLAND 

THRESHOLDED 

BOTTOMLAND 

219 

9 3. ( 

204 

1 

0 

0 

14 

OPEN RANGE 

217 

98.0 

0 

213 

0 

1 

3 

WATER 

399 

98.0 

0 

0 

391 

0 

8 

WOODLAND 

568 

98.0 

1 

4 

0 

556 

7 


Overall Classification Accuracy = 96.8 percent 


BOTTOMLAND : Bottomland 


OPEN RANGE; 


Bermuda, class A (C) ; Bermuda, class E (D) ; Abandoned cropland, 
class A (W) ; Abandoned cropland, class B (X) 


WOODLAND 

WATER 


Mesquite (M) ; Bottomland hardwood (K) ; Dense oak (0) ; Sparse 
oak (P) 

Water (Z) 


I 
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WOODLAND Mesquite (M) i Bottomland hardwood (H) t Dense oak (0) i Sparse oak (P) 
WATER 1 Water (Z) 














































TABLE I -XXI I I.- SNOOK- 1/2 FIRST-ORDER CLASSIFICATION PERFORMANCE 

SUMMARY OF TEST CATEGORIES 


CATEGORY 

NO. OF 
SAMPLES 

% 

CORRECT 

DISTRIBUTION OF SAMPLES 

BOTTOMLAND 

OPEN RANGE 

WATER 

WOODLAND 

THRESHOLDED 

OPEN RANGE 

245 

39.9 

1 

95 

0 

36 

113 

WATER 

119 

98.0 

0 

0 

117 

0 

2 

WOODLAND 

247 

77.0 

0 

0 

0 

191 

56 


Overall Classification Accuracy = 71.3 percent 


BOTTOMLAND: Bottomland (B) 

OPEN RANGE: Bermuda, class A (C) ; Bermuda, class B (D) ; Abandoned cropland, 

class A (W) ; Abandoned cropland, class B (X) 


WOODLAND 

WATER 


Mesquite (M) ; Bottomland hardwood (H) ; Dense oak (O) ; Sparse 
oak (P) 

Water (Z) 


-29 
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TABLE I-XXV.- SNOOK- 1/2 FIRST-ORDER CLASSIFICATION PERFORMANCE 

SUMMARY OF TRAINING FIELDS 


CLASS 

NO. OF 
SAMPLES 

% 

CORRECT 

DISTRIBUTION OF SAMPLES 

BOTTOMLAND 

OPEN RANGE 


WOObUND 

THREShOLDEb 

B 

219 

93.0 

204 

1 

0 

0 

14 

c 

56 

100.00 

0 

56 

0 

0 

0 

D 

94 

98.0 

n 

92 

0 

0 

2 

H 

170 

94.0 

1 

1 

0 

166 

2 

M 

93 

95.0 

0 

3 

0 

88 

2 

0 

250 

99.0 

0 

0 

0 

248 

2 

P 

55 

98.0 

0 

0 

0 

54 

1 

w 

42 

95.0 

0 

40 

0 

1 

1 

X 

25 

100.0 

0 

25 

0 

0 

0 

z 

399 

98.0 

0 

0 

391 

0 

8 


Overall Classification A'icuracy “ 97.0 percent 


BOTTOMUU4D: 
OPEN RANGE: 

WOODLAND 1 


Bottonland IB) 

Bermuda, class A (Or Bermuda, class B (D) r Abandoned cropland, class A (W) t Abandoned 
cropland, class B (XI 

Mesquite (M) t Bottomland hatdwood (H) t Dense oak (0) r Sparse oak (P) 

Water (Z) 


WATER 







I-or.super vised classification pap of Sno 
:if| tiic dSOCLS clustering program (10 major 
; r-c f'ar.eter "valac-s im tais ISOCLS raa v/erc 
1 . 2 , aTDM\>; “ 3.0, f§TilP; ;s;: IS). leference 
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TABLE J-I.- SNOCK-1 DATA INTERPRETATION AND STATISTICS 
OF THE CLUSTERS IN FIGURE J-1 


J-2 




































































lervised classi 
OCLS clusterin^ 
in this ISOCLS 
CLS = 5, ISTOP 


Non 

sup 

using the 
leter valu 
.X * 4.5, 

IS^ 

es 

MAX* 
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TABLE -7-1 1.- SNOOK- 1 DATA INTERPRETATION AND 
STATISTICS OF THE CLUSTERS IN FIGURE J-2 


Cluster 

Color 

Red 

Dark 

Brown 

Bluish 

Green 

Blue 

Yellow 


mean 

mean y^ 

mean y^ 

mean y/ 

mean y/ 

Statistics 

dev. 

/St. dev. 

yst. dev. 

yst . dev . 

y/st . dev . 


33.7 

30.4 y/ 

34.7 y/ 

24.3 y/ 

35.1 y/ 

Channel 1 

9.0 

y/ 1.8 

2.3 

/ 9.3 

/ 

y^ 4.0 


mm 

/ 

22. A y/ 


15.4 y/ 

29.0 y/ 

Channel 2 


y/ 2.4 

u 

y 6.2 

yy/ 7.0 


69.0 

37.8 

41.3 y^ 

10.0 

49.9 y^ 

Channel 3 

9.3 

yy/ 4.2 

y^ 4.6 

y/ 6.5 

y/ 3.5 


38.2 y/ 

20.9 y/ 

21.9 y^ 

3.9 y^ 

26.9 y^ 

Channel 4 

6.0 

y/ 2.3 

y/ 2.6 

// 3.3 

y/ 2.3 

Population 

6,048 

49,970 

38,829 

2,878 

24,775 

Interpre- 

tation 

Cropland , 
very 
vigor- 
ously 
growing 

Woods , 
post oak 

Non- 
wooded 
open 
range ; 
grass- 
land 

Water 

(Lake 

Somer- 

ville) 

Nonwooded 

open 

range and 
vegetated 
land , 
brighter 
in IR 
than the 
green 
cluster 









































ure J-3.“ Honsupervised cl 
sing the ISOCLS clustering 
arameter values in this IS 
TDMAX = 4 . 5 , flA.XCLS = 8,1 



TABLE J-III.- SNOOK-1 DATA INTERPRETATION AND STATISTICS 
OF THE CLUSTERS IN FIGURE J-3 



C| 
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6 \ 






















































Figure J-4,- Classification map of Snook~l data _ obtained by 
a manipulation of the clustering results of figure J-2 . 
[Spectrally similar clusters are grouped together to form 
(1) vigorously growing cropland (red ) , (2) water (blue) , 
(3) other , including rangeland (off-pink) ], 



Figure J-5.- First-order classification map of Snook-1 data 
obtained by manipulation ox the clustering results of 
figure J-3 [clusters that belong to the same classes are 
grouped together to form (1) vigorously growing cropland 
(red ) , (2) water (blue) , (3j clouds (white ) , (4) rangelan 
woodland (brown) , (5) rangeland, open rangeland and grass 

land (green) ] . 
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gure J-6.- Monsupervised classification map 
using the ISOCLS clustering program (10 clus 
values in this ISOCLS run v;ere DLMIN = 0 , ST 
KAKCLS » 10, ISTOP =10). Reference Table J 



'ABLE J-IV.- SNOOK-2 DATA INTERPRETATION AND STATISTICS OF 

THE CLUSTERS IN FIGURE J-6 


J-10 
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Figure J-7.- Monsupervised class! f icaticn map of SB-1 data 

using the ISOCLS clustering program (eight major cluster 
parameter values in this ISOCLS run were DLMIN =0, 
STK-ax * 4.5, HAXCLS * 15, ISTOP = 10). Reference Table 


J-12 
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TABLE J-VI.- SB-2 DATA INTERPRETATION AND STATISTICS OF 
THE CLUSTERS IN FIGURE J-0 
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